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ABSTRACT 


In  May  1964,  a  USAAVLABS-sponsored  turbine/material  research 
program  was  initiated  by  the  AiResearch  Manufacturing  Company 
of  Arizona.  The  overall  program  objective  was  to  perform 
design  studies  and  materials  investigations  to  advance  tech¬ 
nology  for  a  2400°F  uncooled  turbine  component.  The  effort 
was  to  culminate  in  preliminary  design  of  an  uncooled  turbine 
stage . 

The  design  activities  of  the  uncooled-turbine  program  reported 
in  this  document  (Volume  II  of  the  three-volume  final  report) 
were  conducted  concurrently  with  a  materials  research  activity 
to  document  the  material-property  values  of  IN-100  and 
AiResist  13  and  also  to  develop  a  high-temperature  turbine 
material  from  intermetallic  compositions  of  beryllide.  The 
mechanical  analyses  conducted  for  the  turbine  design  were 
based  on  estimates  of  tne  IN-100  and  beryllide  material  prop¬ 
erties. 

Preliminary  design  studies  were  conducted  to  define  a  number 
of  feasible  aerodynamic  turbine  designs  for  efficiency  com¬ 
parisons  and  centrifugal  blade  stress  analyses.  From  these 
studies,  a  candidate  turbine  design  was  selected.  Preliminary 
disk-shape  studies,  thermal  analyses,  and  centrifugal  and 
thermal  elastic  stress  analyses  were  conducted  for  each  rotor 
disk  of  the  candidate  turbine. 

The  design  determined  for  the  uncooled  turbine  evaluations  was 
a  three-stage  axial-flow  turbine  with  an  equal  work  split 
among  the  stages.  Each  turbine  rotor  was  designed  for  40 
blades,  with  the  first  2  stages  using  hollow  blades  of 
beryllide  material  attached  to  IN-100  rotor  disks  by  means  of 
a  2- lobe  fir-tree  attachment  method.  The  third-stage  turbine 
rotor  was  designed  to  be  a  one-piece  casting  of  IN- 100  with 
integral  blades.  The  first-stage  stator  was  designed  as  a 
three-piece-per-vane  configuration  consisting  of  a  platform, 
hollow  vane,  and  outer  shroud— each  piece  to  be  fabricated 
from  beryllide  material.  It  was  anticipated  that  the  second- 
stage  stator  design  would  be  based  on  the  use  of  beryllide 
material  and  that  the  third-stage  stator  material  would  be 
IN- 100.  However,  design  studies  were  not  conducted  for  either 
the  second-  or  third-stage  stators.  Additional  design  fea¬ 
tures  include  a  turbine  inlet  temperature  (TIT)  of  2400 °F,  a 
rotational  speed  of  38,000  rpm,  an  airflow  rate  of  5.0  pounds 
per  second,  a  rotor  disk  cooling-air  flow  rate  of  2.0  percent, 
a  turbine  pressure  ratio  of  9.5:1,  and  a  predicted  totai-to- 
total  turbine  efficiency  of  88  percent. 
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The  feasibility  of  the  uncooied  turbine  design  was  evaluated 
by  mechanical  and  thermal  studies  made  on  the  first-stage 
turbine  components.  Extensive  thermal  and  stress  studies 
were  conducted  to  evaluate  the  mechanical  configurations  of 
ths  rotor  blade,,  disk/blade  attachment,  and  disk,  as  well  as 
the  stator  vane  and  vane  attachment. 

Design  studies  for  the  uncooled  turbine  were  terminated  when 
it  became  apra cent  that  the  proposed  beryl lide  turbine  mate¬ 
rial  was  inadequate  to  satisfy  design  requirements.  As  re¬ 
ported  In  Volume  I  (Materials  Research)  of  the  final  report, 
the  beryllide  compositions  that  were  developed  and  tested 
lacked  sufficient  stress-rupture  strength  for  rotor  applica¬ 
tion  and  did  not  possess  adequate  ductility  for  either  the 
stator  or  the  rotor  application. 


FOREWORD 


The  effort  reported  herein  was  conducted  by  the  AiResearch 
Manufacturing  Company  under  the  terms  of  Contract 
DA  44-177-AMC-183  (T) ,  Task  1G162204A01409. 

The  analytical,  experimental,  and  development  efforts  that 
were  conducted  under  the  uncooled  turbine  program  are  describ¬ 
ed  in  the  three  volumes  of  the  final  report. 

Volume  I  presents  the  results  of  material  investigations 
(both  USAAVLABS-  and  company-sponsored)  that  were  conducted 
to  develop  a  material  that  would  be  suitable  for  the  turbine 
components  of  a  2400 °F  uncooled  turbine. 

Volume  II  presents  the  results  of  the  aerodynamic,  thermo¬ 
dynamic,  and  mechanical  design  activities  that  were  conducted 
for  the  design  of  a  2400°F  uncooled  turbine. 

Volume  III  presents  the  results  of  a  test-evaluation  program 
that  wrs  conducted  to  determine  the  feasibility  of  a  fluidic 
temperature- sensing  system  for  the  measurement  of  TIT*s  in  a 
2400°F  gas  turbine  engine. 

The  manager  of  the  Small  Gas  Turbine  Engine  Project  was  Mr. 

D.  G.  Furst,  and  the  program  manager  of  this  uncooled-turbine 
program  was  Mr.  F.  Weber.  Principal  investigators  were 
Messrs.  R.  W.  Elliott,  J.  L.  Helrobrecht,  G.  E.  Lindstrom, 

J.  J.  Rebcske,  and  A.  S.  Wilson.  The  principal  engineering 
consultants  were  Messrs.  R.  O.  Bullock  and  E.  L.  Wheeler.  The 
overall  guidance  and  technical  direction  provided  by  Messrs. 

J.  White,  H.  Morrow,  L.  Bell,  and  D.  Cale  of  USAAVLABS  are 
gratefully  acknowledged. 
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INTRODUCTION 


This  document  describes  the  aerodynamic,  thermodynamic,  and 
mechanical  design  activities  that  were  conducted  to  evolve  an 
acceptable  unccoled- turbine  design  with  the  following  charac¬ 
teristics: 

1.  TIT  -  2400 °F 

2.  Turbine  airflow  -  approximately  5.0  pounds 
per  sec 

3.  Turbine  efficiency  total  to  total  (T-T)  -  38  percent, 
total  to  static  (T-S)  -  85  percent 

4.  Turbine  stages  -  3  to  5 

5.  Turbine  rotational  speed  -  38,000  to  46,000  rpm 

6.  Turbine  configuration  -  axial-flow  rotors  on 
a  common  shaft 

7.  Turbine  component  materials  -  IN-100,  AiResist  13 
and  an  intermetallic  beryliide  composition  that 
was  concurrently  being  developed 

It  was  expected  that  the  turbine  would  eventually  be  used  in 
a  small  gas  turbine  driving  an  8:1  to  10:1  compressor  and 
furnishing  power  for  single-shaft  engine  application.  The  ob¬ 
jectives  of  this  activity  were  to  establish  the  best  design 
consistent  with  high  efficiency,  mechanical  considerations, 
and  practicality  of  manufacture  for  the  uncooled  turbine. 

Table  I  presents  a  sequential  logic  flow  chart  that  portrays 
the  scope  of  the  overall  design  activities,  shows  the  inter¬ 
relationship  of  the  individual  tasks,  and,  in  addition,  lists 
the  numbers  of  paragraphs  in  which  applicable  discussion  can 
be  found  in  this  document. 
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PRELIMINARY  DESIGN  STUDIES 


2-1  GENERAL  INFORMATION 

The  objective  of  developing  an  uncooled  turbine  with  design 
features  of  2400°F  TIT,  10:1  pressure  ratio,  and  5.0-pound- 
per- second  airflow  required  a  turbine  design  that  was  opti¬ 
mized  from  a  stress  standpoint  in  addition  to  the  development 
of  an  advanced  turbine  material.  As  reported  in  Volume  I,  a 
materials  research  program  was  being  concurrently  conducted 
with  the  turbine  design  studies  to  develop  a  beryllide  com¬ 
position  with  sufficient  material  properties  for  use  as  a 
high-temperature  turbine  material.  Since  material-property 
design  data  were  not  yet  available  for  the  beryllide  material, 
the  preliminary  design  studies  placed  emphasis  on  reducing 
the  inherent  stress  levels  in  the  turbine  design.  The  result¬ 
ing  stress  levels  then  became  the  objective  of  the  materials 
research  program.  It  was  assumed  that  the  beryllide  composi¬ 
tion  material  when  developed  would  be  adequate  for  use  in  the 
rotor  blades  of  the  first  two  turbine  stages  and  also  for  use 
in  the  first-  and  second-stage  stators. 

Designs  for  the  first-  and  second-stage  rotor  disks,  all 
subsequent  stage  rotors,  and  all  stage  stators  except  for  the 
first  two  stages  were  based  on  the  premise  that  the  turbine 
material  would  be  IN- 100.  Extensive  literature  research  was 
conducted  to  determine  the  present  availability  of  material 
properties  for  IN-100.  The  results  revealed  that  current 
stress-rupture  data  for  IN-100  showed  considerable  scatter. 
The  IN-100  data  scatter  was  found  to  lie  both  above  and  below 
the  current  minimum  stress-to-rupture  values  for  as-cast 
INCO  713C.  Therefore,  it  was  decided  that  INCO  713C  data 
would  be  used  for  the  design  studies  until  such  time  as  reli¬ 
able  IN-100  data  became  available  from  the  IN- 100  materials 
test  program,  reported  in  Volume  I.  Thus,  if  IN- 100  failed 
to  meet  expectations,  it  would  not  be  unreasonable,  from  a 
stress  consideration,  to  substitute  as-cast  INCO  713C. 

The  preliminary  design  studies  were  directed  to  define  a  group 
of  candidate  aerodynamic  turbine  designs,  from  which  one 
design  could  be  selected  on  the  basis  of  centrifugal  blade 
stress  evaluations  and  attainable  turbine  efficiency  compari¬ 
sons. 
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It  was  believed  that  the  last  stage  of  the  turbine  would  have 
the  roost  critical  centrifugal  blade  stress  levels.  This 
opinion  was  based  on  (1)  the  assumption  that  the  conventional 
material  of  the  last-stage  components  would  be  more  vulnerable 
than  the  beryl lide  material  of  the  upstream  stages  with  re¬ 
spect  to  estimated  stress  conditions,  and  (2) ,  the  fact  that 
the  last-stage  annular  area  requirement  results  in  greater 
blade  length*  Therefore,  the  preliminary  centrifugal  blade 
stress  evaluations  were  conducted  on  the  exit  stage  only. 

To  initiate  the  preliminary  design  studies,  cycle  calculations 
were  performed  to  establish  the  turbine  exit  temperature, 
average  gas  properties,  total-to-ambient  pressure  ratio,  and 
turbine  inlet  equivalent  flow.  The  last-stage  exit  critical 
velocity  and  diffuser  static-pressure-recovery  coefficient 
were  then  selected,  with  consideration  given  to  turbine  effi¬ 
ciency  loss  and  to  centrifugal  stresses  related  to  the  result¬ 
ing  exit  annulus  area  established.  The  number  of  turbine 
stages,  stage  work  split,  and  rotational  rpm  were  then  varied 
to  establish  the  physical  size  of  the  turbines  and  the  at¬ 
tendant  vector  diagrams  and  hub  speeds  (i.e.,  rotor  disk  rim 
speeds) .  Nine  candidate  preliminary  aerodynamic  turbine 
designs  were  defined  by  these  procedures  for  further 
evaluation,  comparisons,  and  finally,  selection  of  one  design 
to  critique  for  final  design  studies. 

2.2  CYCLE  CALCULATIONS 

Cycle  calculations  were  made  to  establish  the  overall  require¬ 
ments  of  the  turbine  in  terms  of  turbine  flow  and  work  out¬ 
put.  A  single-spool  engine  with  the  following  operational 
cycle  for  a  sea-level  standard  day  was  selected: 

1.  Compressor  pressure  ratio,  10:1 

2.  Compressor  airflow,  5.0  pounds  per  second 

3.  Compressor  efficiency,  82  percent  (T-T) 

4.  Turbine  rotor  disk  cooling  airflow,  2  percent  of 
compressor  airflow 

5.  Burner  pressure  loss  Ap/p,  0.05 

6.  Burner  efficiency,  98  percent 

7.  Lower  heating  value,  18,400  Btu  per  pound 

8.  TIT,  2400 °F 
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9.  Turbine  efficiency,  -  (T-T)  88  percent, 

(T-S)  85  percent 

With  use  of  the  above-listed  cycle  parameters,  the  overall 
turbine  design  requirements  were  calculated  to  be: 

1.  Turbine  exit  temperature,  1381°F 

2.  Turbine  inlet  equivalent  flow,  1.247  pounds 
per  second 

3.  Turbine  pressure  ratio  (total- to-ambient) ,  9.5:1 

In  addition,  a  speed  range  of  30,000  to  42,000  rpm  was  speci¬ 
fied.  It  was  judged  that  feasible  compressor  designs  could 
be  obtained  in  this  speed  range,  with  the  higher  speeds  being 
more  desirable  from  the  compressor  design  point  of  view  and 
the  lower  speeds  being  more  desirable  from  a  turbine  stress 
point  of  view. 

With  these  general  turbine  requirements  defined  and  the  speed 
range  specified,  turbine  studies  were  made  to  determine  the 
quantitative  dependence  of  turbine  efficiency  and  last-stage 
centrifugal  blade  stress  on  the  choice  of  exit  axial  critical 
velocity  ratio,  exhaust  diffuser  static-pressure-recovery 
coefficient,  and  rotative  speed. 

2.3  LAST-STAGE  EXIT  CRITICAL  VELOCITY  RATIO 

For  the  subject  uncooled- turbine  design,  selection  of  the  exit- 
stage  critical  velocity  ratio  was  made  by  determining  the  best 
compromise  between  exhaust  loss  in  terras  of  attainable  turbine 
efficiency  and  the  centrifugal  stresses  directly  related  to 
the  rotor  exit  annulus  area.  The  following  subparagraphs 
detail  the  studies -  that  were  made  to  select  an  exit-stage 
velocity  ratio  of  0.45. 

2.3.1  Effect  of  Last-Stacfe  Exit  Critical  Velocity  Ratio 
on  Attainable  Turbine  Efficiency 

Any  kinetic  energy  leaving  the  exhaust  plane  of  a  shaft-power 
engine  diffuser  exit  is  considered  to  be  a  loss.  This  loss 
can  be  expressed  as  the  ratio  of  two  isentropic  available 
energies  across  the  turbine.  First,  if  it  is  assumed  that  a 
perfect  diffuser  and  infinite  area  are  present,  the  velocity 
at  the  exit  plane  of  the  engine  can  be  made  zero.  For  this 
case,  the  available  energy  to  the  turbine  is  based  on  the 
inlet  total- to- ambient  pressure.  The  second  is  the  case  where 
a  real  diffuser  is  considered.  The  isentropic  available 
energy  here  will  be  defined  by  the  inlet  total-to-exit 
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total  pressure  ratio  and  is  smaller  than  the  first  case  by  the 
real  diffuser  losses.  This  pressure  ratio  can  in  turn  be  re¬ 
lated  to  the  axial  critical  velocit'*  ratio  at  the  exit  of  the 
turbine  and  a  diffuser  static-pressure-recovery  coefficient. 

For  a  given  TIT  and  average  gas  properties,  the  ratio  of  the 
isentropic  available  energies  can  be  expressed  as 


*** (total -to tal) 

^ (total-ambient) 


(1) 


The  decrease  in  attainable  efficiency  on  an  isentropic  basis 
is  therefore 


^ (total  to  total) 
^ (total  to  static) 


(2) 


For  a  fixed  overall  pressure  ratio  — —  ,  the  total-to-total 
pressure  ratio  is  Famb 


Pexit  _  ( Pamb\  /Pexit\ 
p'in  \pin  /  \pamb  / 


(3) 


The  diffuser  static  pressure  recovery  as  the  ratio  of  the 
static  pressure  rise  to  the  difference  in  total-to-static 
pressure  at  the  turbine  exit  is  defined  as 


P  —  P 
aiab  exit 

"iP  I  p 

exit  exit 


(4) 


6 


and 


p;Xit 

Pamb 


1  + 


(P'/P) 


exit 


(5) 


For  the  condition  of  no  exit  swirl,  the  ratio  of  the  total-to 
static  pressure  at  the  turbine  rotor  exit  is  a  function  of 
only  the  exit  axial  critical  velocity  ratio. 


(6) 


Figure  l  shows  the  decrease  in  attainable  efficiency  as  a 
function  of  exit  axial  critical  velocity  ratio  for  an  overall 
total- to-anibient  pressure  ratio  of  9.5  and  a  diffuser  static- 
pressure-recovery  coefficient  Rp  =  0.6.  Experience  has  shown 
that  this  value  of  Rp  is  about  the  maximum  value  that  can  be 
expected  for  an  annular  diffuser  operating  at  the  turbine 
exit. 

2.3.2  Effect  of  Last-Stage  Exit  Axial  Velocity  and  Rotative 
Speed  on  Rotor  Blade  centrifugal  Stress 

The  average  blade-root  centrifugal  stress  may  be  represented 
as  a  function  of  the  rotative  speed  and  exit  annular  area  by 
the  following  equation  for  centrifugal  loading: 


„  -h'^X 

“  24nq 


(7) 


where  PM  =  blade  material  density,  Ib/cu  in. 

$  =  blade  taper  factor,  dimensionless 
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Figure  1.  Loss  in  Attainable  Efficiency. 
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tu  =  rotational  velocity,  radians/sec 

Aa  -  rotor  annulus  area,  sq  in. 

g  -  acceleration  of  gravity,  ft/sec 

o  =  average  blade  root  stress,  psi 

The  annular  area  in  this  equation  can  be  expressed  (again  for 
the  case  of  no  exit  swirl  component  of  velocity)  as  a  func¬ 
tion  of  the  rotor  exit  axial  critical  velocity  ratio,  the 
diffuser  static-pressure-recovery  coefficient,  and  the  equiv¬ 
alent  flow  required  at  the  last-stage  rotor  exit.  The  exit 
equivalent  flow  required  is  expressed  as 


w/T 

~~i —  (8) 

re 


where  w 
6re 
e 


actual  exit  flow,  Ib/sec 
Pre/Parob  ‘  f(^)  '  and  *! 


T’/T  . 
e  amb 


The  equivalent  flow  per  unit  of  area  is  expressed  in  terras  of 


—  as  follows : 
cr 


w/J 


where 


Y+X  \  *'CJ  reJ  V'J 


1 


acr.  sta  -  im  9  itTstd  j 


cr  re 

1/2 


p*  std  a*  std 
cr 


(9) 


std 


°  std  R  T 


std 
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For  the  assumed  static-pressure-recovery  coefficient  of  0.6, 

w/jr  v 

Figure  2  shows  6 re  and  as  a  function  of  The  required 

A  cr 

annular  area  then  is  simply 


W-/TT  r-jr- 

a  _ _ 5  /  yZdJL- 

aa  6  7  a7  y 

re  a 


(10) 


For  assumed  values  of  rpm,  blade  material  density,  and  blade 
stress  taper  factor,  it  is  now  possible  to  obtain  the  blade 

root  stresses  as  a  function  of  ^F—  at  the  rotor  exit.  Such  a 

acr  / 

plot  for  blade  root  stress  is  shown  .in  Figure  3  for  a  range 
of  rpm,  a  conventional  blade  material  density  of  0.3  pound 
per  cubic  inch  (the  density  of  ZNCO  713C),  a  diffuser  recovery 
factor  of  0.6,  and  an  assumed  blade  taper  factor  of  0.73. 

Figure  3  shows  that  for  any  rotative  speed,  the  rotor  blade 
centrifugal  stress  can  be  significantly  reduced  by  going  to 
a  higher  rotor  exit  critical  velocity  ratio.  Figure  1,  how¬ 
ever,  shows  that  there  is  an  accompanying  decrease  in  attain¬ 
able  turbine  efficiency  associated  with  increasing  exit 
velocity.  For  this  cycle,  a  two- point  decrease  in  turbine 
efficiency  would  decrease  the  specific  power  by  8.5  percent 
or,  for  a  given  specific  power,  increase  the  inlet  tempera¬ 
ture  by  60®F.  Based  on  these  considerations,  a  value  of  0.45 
was  selected  as  a  reasonable  compromise  for  the  rotor  exit 
axial  critical  velocity  ratio  despite  the  indicated  higher 
blade  stresses  in  Figure  3. 

It  would  also  appear  from  Figure  3  that  the  turbine  speed  is 

V„ 

uniquely  defined  by  the  specification  of  T—  and  a  blade 

acr 

centrifugal  stress,  and  to  a  large  extent  this  is  true.  How¬ 
ever,  the  blade  taper  factor  ♦  assumes  a  linear  blade  profile 
area  variation  with  respect  to  radius  and  can  have  an  effect 
on  the  calculated  centrifugal  stress.  For  instance,  a  loga¬ 
rithmic  or  hyperbolic  blade  taper  would  result  in  .reduced 
centrifugal  stress  for  the  same  exit  velocity  and  rpm.  For 
this  reason,  more  exact  knowledge  of  the  turbine  blading  is 
required  before  the  centrifugal  stresses  can  be  known  for  a 
given  exit  velocity  and  rpra. 
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2.3.3  Effect  on  Disk  Stress  of  Stage  work  Coefficient  and 
Number  of  Stages 

The  average  tangential  disk  stress  is  dependent  upon  (among 
other  factors)  the  rim  load,  which  is  the  product  of  the 
blade  centrifugal  stress  and  the  total  blade  profile  area  at 
the  rim,  and  the  rim  velocity  or  turbine  hub  speed.  There¬ 
fore,  lower  hub  speed  would  result  in  a  lower  disk  stress. 
The  hub  speed  of  the  turbine  for  a  given  work  requirement  is 
dependent  upon  the  product  of  the  number  of  stages  and  the 
stage  work  coefficient.  For  a  turbine  of  free  vortex  design 
and  constant  hub  radius,  the  total  work  output  can  be  ex¬ 
pressed  by: 


n 
a 

Required  work  =  gJAK  =  u 

n 

i 


(id 


where 


or 


Xg  =  AVu/ur  =  stage  work  coefficient 

Av  =  change  in  tangential  component 
of  velocity  through  the  rotor 

U„  =  hub  tangential  velocity,  ft 
H  per  sec 

n  --  number  of  stages 


n 

required  work/ 

1 


X 


s 


The  work  requirement  can  be  met  at  a  given  hub  velocity  either 
by  increasing  the  stage  work  coefficient  or  by  increasing  the 
number  of  stages  with  a  given  work  coefficient. 

Quantitative  evaluations  of  these  effects  on  disk  stress, 
however,  required  more  detailed  knowledge  of  the  turbine 
geometry  for  a  given  number  of  stages  and  individual  stage 
work  coefficients. 
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2.4  CHOICE  OP  STAGE  WORK  COEFFICIENTS 
- . —  — . ■— - -  ■  1  ■  . . —  1  ■■  ■  ■ 

v 

Figure  4  shows  a  typical  vector  diagram  (nondimens iona 1 i zed 
by  the  hub  speed)  and  station  identification.  For  a  constant 
hub  radius  from  Station  2  to  3  and  with  assumption  of  a  free- 
vortex  distribution  of  velocity  with  respect  to  radius,  the 
stage  work  coefficient  can  be  expressed  (from  Euler's  work 
equation)  as 


Xstage  X2,H  “  *3,H 


where 


(13) 


It  is  also  apparent  from  the  vector  diagram  that 


(14) 


and  therefore 


+  1 


(15) 


and 


(16) 
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Figure  4.  Typical  Hub  Vector  Diagran 
and  Station  Identification. 
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If,  for  the  sake  of  simplicity,  constant  axial  velocity  is 
assumed  at  Stations  2  and  3,  then  the  ratio  of  the  relative 
tangential  velocity  is  the  same  as  the  ratio  of  the  relative 
velocity.  Letting  K  be  this  ratio,  we  find  that 


(negative  sign  required 
for  change  in  vector  (17) 
direction) 


where 


For  K  >  i,  there  is  an  accelerating  flow  through  the  rotor — 
i.e.,  a  reaction  design.  For  K  -  1.0,  the  relative  velocities 
entering  or  leaving  the  rotor  are  equal — i.e.,  an  impulse 
design.  For  K  <  1.0,  there  is  a  velocity  deceleration  across 
the  rotor — i.e.,  a  negative  reaction  design,  without  proper 
boundary-layer  control,  velocity  deceleration  across  the 
rotor  is  to  be  avoided,  since  separation  is  likely  to  occur, 
with  large  attendant  losses.  In  terms  of  the  relative  tan¬ 
gential  velocities,  then,  the  stage  work  coefficient  is 


'  (^h  a  J  K>  <18) 


and 


1  f  1  +  K 


Now,  if  we  specify  zero  exit  swirl  (13  =0),  then  \s  =  \2,h 
and  rotor  hub  impulse  conditions  (i.e.,  K  =  1.0),  it  then 
follows  that 


2.0 


(20) 
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For  Xs  >2.0  and  K  =  1.0,  X3  takes  on  negative  values.  Which 
decreases  the  reaction  available  and  increases  the  turning 
requirement  on  the  following  downstream  stator.  If  the  stage 
is  a  last  stage,  the  kinetic  energy  represented  by  the  tan¬ 
gential  component  of  velocity  cannot  be  recovered  in  an 
annular  diffuser  and  is  therefore  a  loss.  For  these  reasons, 
turbines  with  stage  work  coefficients  greater  than  2.0 
generally  exhibit  lower  aerodynamic  efficiencies  than  those 
with  stage  work  coefficients  equal  to,  or  less  than,  2.0. 

Additionally,  equal  entering  and  leaving  relative  velocity  at 
the  rotor  hub  {K  =  1.0 — impulse  design)  is  usually  accepted 
as  being  a  limiting  aerodynamic  condition  that  is  consistent 
with  good  efficiency.  Impulse  designs  can  also  produce 
higher  work  output  for  a  given  rotor  blade  hub  speed  than 
cam  reaction  designs. 

2.5  CANDIDATE  TURBINE  DESIGNS  FOR  STRESS  AND 
EFFICIENCY  EVALUATIONS 

Nine  candidate  turbine  designs  were  prepared  for  centrifugal 
blade  stress  and  efficiency  analyses,  from  which  the  most  ad¬ 
vantageous  design  was  to  be  selected  for  final  design  studies. 
Based  on  the  reasoning  presented  in  the  preceding  paragraphs, 
each  of  the  designs  had  an  exit-stage  critical  velocity  ratio 
of  0.45  and  impulse  conditions  specified  at  the  rotor  blade 
hubs.  Additionally,  most  of  the  designs  had  stage  work  coef¬ 
ficients  of  2.0. 

Table  II  presents  a  summary  of  the  candidate  turbine  designs 
whose  gas  paths  were  defined  in  detail  and  shows  the  speed 
range,  number  of  stages,  and  stage  work  coefficients  con¬ 
sidered. 

Subsequent  paragraphs  report  the  stress  and  efficiency  analy¬ 
ses  that  were  conducted  to  select  the  most  advantageous 
design. 

2.5.1  Stress  Evaluations  of  the  Candidate  Turbine  Designs 

Before  the  blade  centrifugal  stress,  stress-rupture  life,  and 
average  tangential  stress  in  the  rotor  disk  can  be  accurately 
determined,  a  number  of  variables  must  be  known.  Among  these 
variables  are  an  accurate  and  complete  geometrical  description 
of  the  blade,  blade  metal  temperatures,  material  properties, 
and  the  method  by  which  the  blade  is  attached  to  the  disk. 
Therefore,  prior  to  initiation  of  stress  evaluations  of  the 
candidate  designs,  efforts  were  conducted  to  either  determine 
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TABLE  II,  NINE  CANDIDATE  PRELIMINARY 
AERODYNAMIC  DESIGNS 


1 

2 

Stage 

3 

Number 

4 

5 

1 

| 

Stage  work 

Coefficients 

a« 

2.0 

2.0 

2.0 

— 

- 

6.0 

►a 

• 

o 

2.0 

2.0 

2.0 

w*» 

8.0 

2.0 

2.0 

2.0 

2.0 

2.0 

io.  e 

2.0 

2.0 

2.0 

- 

- 

6.0 

2.0 

2.0 

2.0 

2.0 

- 

8.0 

2.0 

2.0 

2.0 

- 

•»* 

6.0 

2.0 

! 

2.0 

2.0 

2.0 

- 

8.0 

i  3.0 

3.0 

2.0 

mm 

mm 

8.0 

3.0 

3.0 

2.0 

- 

- 

8.0 

Turbine  Design 
Identification 
Number 


A-338-5 

A-438-5 

A-538~5 

A-342-5 

A-442-4 

A-346-4 

A-446-5 

B-338-9 

B-338-10* 


38,000 

38,000 

38,000 

42,000 

42,000 

46,000 

46,000 

38,000 

38,000 


*Same  blade  heights  as  Turbina  A-338-5. 
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or  estimate  the  required  information  for  the  subject  var¬ 
iables.  The  following  subparagraphs  report  these  efforts  and 
also  the  results  of  the  stress  evaluations. 

2.5. 1.1  Blade  Geometry 

A  sizable  man-hour  effort  is  required  to  define  a  complete 
and  suitable  three-dimensional  blade  shape.  The  procedure 
is  to  specify  a  two-dimensional  profile  shape  of  a  given 
radius#  determine  the  surface  velocity  distribution  from  a 
channel  flow  solution,  and  subject  the  profile  shape  to 
perturbation  until  a  suitable  velocity  distribution  is  ob¬ 
tained.  This  procedure  is  repeated  at  a  number  of  specified 
radii  between  the  hub  and  tip  shrouds  of  the  flow  passage. 
The  two-dimensional  profiles  are  then  stacked  with  their 
center  of  gravity  on  a  radial  line  and  faired  radially  from 
hub  to  tip.  If  a  smooth  three-dimensional  fairing  does  not 
result,  the  procedure  is  repeated  until  a  smooth  three- 
dimensional  shape  is  evolved.  This  blade  is  then  further 
analyzed  with  a  complete  three-dimensional  flow  solution. 

Thus,  the  work  involved  in  completing  a  three-dimensional 
blade  shape  for  each  of  the  nine  candidate  turbine  designs 
for  stress  analyses  would  have  been  prohibitive  for  this 
study.  Therefore,  a  method  was  conceived  to  generate  two- 
dimensional  blade  profiles  at  each  radial  station  uniquely 
specified  by  the  entering  and  leaving  angles  of  the  local 
vector  diagram  and  Zweifel's*  solidity  requirements.  The 
blade  section  areas  were  then  calculated  and  integrated,  to¬ 
gether  with  their  radial  moments,  to  obtain  blade  volumes 
and  radial  centroids  as  reported  in  Appendix  1  of  this  docu¬ 
ment. 


2. 5*1. 2  Blade  Metal  Temperatures 

In  addition  to  the  centrifugal  blade  stresses  calculated 
based  on  the  blade  geometry  study  discussed  above,  it  was 
necessary  to  determine  the  blade  metal  temperatures  before 
the  blade  stress-rupture  life  could  be  determined.  For  this 
study,  the  blade  metal  temperatures  were  assumed  to  be  equal 


*  ZWaifel,  0.#  OPTIMUM  BLADE  PITCH  FOR  TURBOMACHINES  WITH 
SPECIAL  REFERENCE  TO  BLADES  OF  GREAT  CURVATURE#  Brown 
Bovsri  Review#  Vol.  37#  No.  12,  December  1945,  pp.  436- 
444. 
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to  the  adiabatic  wail  temperature  at  each  radial  section,  ft 
"flat"  temperature  profile  {total  gas  temperature  constant 
across  the  annulus)  was  assumed  at  the  rotor  inlet.  The  adi¬ 
abatic  wall  temperature  (TAW)  was  then  determined  from  the 
static  gas  temperature  (Ts) ,  total  gas  temperature  relative 
to  the  blade  (T») ,  and  an  appropriate  recovery  factor  (R  )  at 
each  radial  section.  The  relationship  between  these  quanti¬ 
ties  is  shown  below: 


W2 


*P*  =  m  +  _ 

R  S  2gJc 


(21) 


W2 


TAW  -  TS  +  Rc  2gJc~  *  Tc  +  (<  ~  T J 


C  R 


(22) 


where  Rc  is  equal  to  the  cube  root  of  the  Prandtl  number  for 
turbulent  flow,  and  the  term  represents  the  kinetic 

energy  of  the  gas  relative  to  the  blade. 

2*5-1-3  Material  Properties  and  Blade  Attachment  Method 

For  the  reasons  given  in  Paragraph  2.1  of  the  preliminary 
design  studies,  the  centrifugal  blade  stress  analyses  were 
conducted  only  on  the  exit  stage  of  the  candidate  turbine 
designs.  The  exit-stage  rotor  designs  for  each  of  the 
candidate  turbines  were  of  a  one-piece  rotor  configuration 
with  integrally  cast  blades. 

INCO  713C  material  property  data  were  used  for  the  analyses. 
Stress-rupture  values  for  INCO  713C  are  presented  in  the 
Larson-Miller  diagram  of  Figure  5.  The  design-life  objective 
of  the  uncooled  turbine  was  200  hours  at  2400 °F  first-stage 
TIT. ^  substituting  a  value  of  200  hours  for  the  time  variable 
in  Figure  5  results  in  the  plot  of  design  stress  level  versus 
temperature  shown  in  Figure  6.  Included  in  Figure  6  is  the 
same  plot  for  IN-100  data  that  was  later  obtained  from  the 
materials  test  program  reported  in  Volume  I. 
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Figure  5.  INCO  713C  Stress-Rupture  valiirs. 
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2.5. 1.4 


Results  of  Centrifugal  Blade  Stress  Evaluations 


•Che  three-stage,  equal-work-split,  38,0QQ-rpra  turbine  (Design 
A-338-5)  was  found  to  have  the  greatest  stress-rupture  life 
(in  excess  of  the  200-hour  design  goal) ;  thus,  from  a  stress 
standpoint,  it  became  the  most  promising  configuration. 
Candidate  turbine  designs  eliminated  from  further  consider¬ 
ation  because  they  did  not  meet  the  design  life  goal  were: 


1.  Design  A-342-5 

2.  Design  A-346-a 

3 .  Design  A-446-5 

Additionally,  it  was  decided  to  eliminate  Design  a -442  -4.  Al¬ 

though  the  last  stage  of  this  four -stage  design  met  the  design 
life  goal,  the  lack  of  adequate  stress  margin  dictated  beryl - 
xide  blading  in  the  first  three  stages.  The  added  complexity 
of  an  additional  stage  with  beryl lide  blading  and  the  addi¬ 
tional  weight  of  a  fourth  stage  were  the  reasons  for  elimi¬ 
nating  this  design  from  further  consideration.  The  remaining 
candidate  turbine  designs  all  seemed  feasible  from  a  stress 
standpoint  and  were  submitted  for  turbine  efficiency  evalua¬ 
tions. 

Representative  data  for  these  centrifugal  blade  stress  analy¬ 
ses  are  shown  in  Figures  7,  8,  and  9  for  Designs  A-338-5, 
A-342-5,  and  A-442-4,  respectively. 

2.5.2  Efficiency  Evaluations  of  the  Candidate  Turbine  Designs 

When  the  candidate  turbine  designs  were  being  defined  for  sub¬ 
mittal  to  stress  analyses,  an  attempt  was  made  to  keep  blade 
heights,  turning  angles,  and  blockage  factors  within  the  re¬ 
gions  of  knswn  good  efficiencies.  However,  it  could  not  be 
expected  that  the  various  designs  would  yield  identical  tur¬ 
bine  efficiencies.  Therefore,  the  candidate  turbine  designs 
that  were  determined  to  be  feasible  from  the  centrifugal  blade 
stress  analyses  (Paragraph  2.5.1)  were  subjected  to  turbine 
efficiency  evaluations  and  comparisons. 

It  was  determined  during  preliminary  valuations  that  the 
five-stage  Turbine  Design  A-538-5  wou-u  not  provide  sufficient 
efficiency  improvements  over  other  candidate  designs  to  off¬ 
set  the  disadvantages  of  the  increased  weight  and  additional 
rotating  hardware.  Therefore,  Turbine  Design  A-538-5  was 
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NOTES 

*.  ROTOR  SPEE0=38,000  PPM,  TIP  SPEED=1630  FPS 
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3.  MATERIAl=INCO  7!3C,  AS  CAST 

4.  INLET  TEMPERATURE=1752  °F 

5.  THIRD  STAGE 

6.  NO.  OF  8LADES=40 


Figured.  Preliminary  Centrifugal  Blade 
Stress  Analysis  for  Grit  Stage 
of  Turbine  Design  A-338-5 . 
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Figure  8,  Preliminary  centrifugal  Blade  Stress 
Analysis  for  Exit  Stage  of  Turbine 
Design  A-342-5. 
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4.  INLET  TEMPERATURE=1672  °F 

5.  FOURTH  STAGE 
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Figure  9.  Preliminary  Centrifugal  Blade  Stress 
Analysis  for  Exit  Stage  of  Turbine 
Design  A-442-4. 
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eliminated  from  further  consideration.  This  elimination,  plus 
eliminations  resulting  from  the  stress  analyses  (Paragraph 
2. 5. 1.4),  left  only  four  of  the  original  nine  designs  as 
candidates : 


1.  Turbine  Design  A-338-5 

2.  Turbine  Design  A-438-5 

3.  Turbine  Design  B-338-9 

4.  Turbina  Design  B-338-10 

preliminary  evaluations  of  the  two  3-nurabered  designs  re¬ 
vealed  that  Design  3-338-10  afforded  greater  efficiency 
potentials.  Therefore,  comparative  efficiency  studies  were 
conducted  on  three  candidate  designs — A-338-5,  A-438-5,  and 
B-338-10.  Appendix  II  presents  a  summary  of  the  methods  that 
were  employed  to  estimate  the  aerodynamic  efficiency  of  the 
turbine  designs. 


The  three  aerodynamic  turbine  designs  that  were  submitted  to 
the  efficiency  studies  incorporated  comparative  solutions  for 
reducing  hub  diameters  and,  consequently,  hub  speeds,  in  order 
to  keep  centrifugal  stresses  of  the  rotors  at  a  minimum.  Each 
of  the  turbines  was  designed  for  a  rotor  speed  of  38,000  rpra. 
Turbine  Design  A-338-5  has  three  impulse  hub  ( \s  =  2.0) 
stages,  with  a  resulting  hub  radius  of  3.319  inches.  Turbine 
Design  A-438-5  has  four  similar  stages,  with  a  resulting  hub 
diameter  of  2.875  inches.  Turbine  Design  B-338-10  attempts 
to  accomplish  four— stage  loading  in  three  stages  by  using 
negative  swirl  in  the  first  two  stages,  and  also  has  a  hub 
diameter  of  2.875  inches. 


Figure  10  shows  comparative  total-to-static  turbine  efficien¬ 
cies  for  these  three  turbines.  It  can  be  seen  that  the  four- 
stage  solution  (A-438-5)  to  lower  hub  radius  results  in  im¬ 
proved  turbine  efficiency,  while  the  heavy  loading  of 
Turbine  Design  B— 338— 10  results  in  lowered  efficiency.  It 
may  be  possible  to  minimize  the  performance  loss  of  the 
highly  loaded  three-stage  turbine  (B-338-10)  by  increasing 
tne  blade  heights  in  the  f  irst  two  stages  (B-338-10  has  th® 
same  blade  heights  as  A-338-5).  However,  the  comparative  ~  . 
-zero-clearance  losses  and  thfe  slight  improvement  In  0.030-  1 

of  A-438~5  over  A-338-5  (due  to  increased 
olade  height)  indicate  chat  Turbine  Design  B-338-10  would 
still  have  inferior  performance  with  blade  heights  as  tail 
as  those  of  A-438-5.  Therefore,  Turbine  Design  B-338-10  and, 

•Consequently,  B-338-9  were  eliminated  from  further  considera¬ 
tion. 
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Figure  10.  Comparison  of  Efficiencies 
for  Three  Candidate  Turbine 
Designs. 


2 . C  SELECTION  OF  PRELIMINARY  TURBINE  DESIGN 

The  two  most  advantageous  candidate  turbine  designs  to  pursue 
for  final  design  activities  Were  A-338-5  and  A-438-5.  In  the 
final  analysis,  it  was  decided  to  use  the  A-338-5  design 
because  of  two  factors:  (1)  the  calculated  centrifugal  blade 
stresses  were  lower  for  Design  A-338-5,  and  (2)  the  belief 
that  the  higher  estimated  efficiency  of  A-438-5  did  not  off¬ 
set.  the  weight  penalty  and  additional  rotational  hardware 
imposed  by  use  of  a  fourth  stage. 

Hie  eight  preliminary  aerodynamic  candidate  turbine  designs 
that  were  determined  to  be  less  advantageous  than  Design 
A-338-5  are  listed  below  with  the  characteristics  that 
eliminated  them  from  further  design  considerations. 


1. 

A-438-5, 

less  stress-rupture  life  than  A-338-5 
and  not  enough  efficiency  improvement 
to  warrant  use  of  a  fourth  stage. 

2. 

A-538-5, 

less  stress-rupture  life  than  A-338-5 
and  not  enough  efficiency  improvement 
to  warrant  use  of  a  fourth  and  fifth 
stage. 

3. 

A-342-5, 

did  not  meet  the  200-nour  stress-rupture 
design  life. 

4. 

A-442-4, 

marginal  200-hour  stress-rupture  design 
life  necessitated  three  turbine  stages 
with  beryllide  blading  plus  the  dis¬ 
advantages  of  four  stages. 

5. 

A- 346-4, 

did  not  meet  the  200-hour  stress-rupture 
design  life. 

6. 

A-446-5, 

did  not  meet  the  200-hour  stress-rupture 
design  life. 

*  7. 

B-338-9, 

less  stress-rupture  life  and  efficiency 
than  A-338-5. 

8. 

B-338-10, 

less  stress-rupture  life  and  efficiency 
than  A-338-5 .  ™ 
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2.7 


PRELIMINARY  DISK  SHAPES 


Two  quantities  remained  to  be  determined  before  the  turbine 
disk  shapes  could  be  defined:  the  ratio  of  mid-blade  thick¬ 
ness  to  axial  chord  (TC/CX)  and  the  number  of  rotor  blades 
per  stage  (see  Appendix  I).  The  subject  quantities  were 
established  by  determining  a  suitable  compromise  between  aero¬ 
dynamic  performance  and  the  centrifugal  stress  levels  ir»  the 
disk,  with  consideration  given  to  prior  turbine  programs  where 
successful  fabrication  and  casting  experience  had  been  gained 
on  similar  integrally  cast  and  nonintegrally  fabricated 
turbine  rotors  of  the  sane  si2e.  The  values  of  TC/CX  that 
were  being  considered  were  0.15  and  0.2G.  The  number  of 
blades  under  consideration  were  30  and  40. 

It  was  determined  that  both  stress  and  efficiency  were  im¬ 
proved  by  the  higher  number  of  blades  (40) .  The  thinner 
blades,  TC/CX.  =  0.15,  were  determined  to  be  more  desirable 
from  a  disk  stress  standpoint  but  resulted  in  a  slight  loss 
of  aerodynamic  performance.  Six  studies  were  then  made  to 
select  the  final  configuration  of  each  disk  (se*=»  Table  III)  - 

For  the  third-stage  disk,  two  configurations  were  studied  and, 
as  shown  in  Table  III,  a  blade  number  of  40  and  a  TC/CX  chord 
ratio  of  0-20  were  selected.  The  thicker  chord  ratio  was 
chosen  in  order  to  gain  the  additional  aerodynamic  performance 
available,  as  prior  experience  with  integrally  cast  rotors 
with  the  same  number  of  blades  and  blade  thickness  and  similar 
operating  environment  had  not  resulted  in  disk  stress  problems. 

As  shown  in  Table  III,  three  configurations  were  studied  for 
the  second- stage  disk  where  both  TC/CX  and  the  number  of 
blades  were  varied  to  determine  the  effect  on  disk  sizing.  To 
assure  blade  spacings  large  enough  to  allow  for  the  antici¬ 
pated  fir-tree  attachment  of  the  blades  to  the  disk  and  to 
attain  the  best  aerodynamic  performance,  the  blade  number 
selected  was  40.  Because  of  stress  considerations  to  the 
attachment,  the  chord  ratio  chosen  was  0.15.  Based  on  the 
results  of  the  studies  performed  on  the  second-stage  disk, 
only  one  first-stage  disk  was  considered.  As  shown  in  Table 
III,  the  blade  number  for  the  first  stage  was  40  and  the  chord 
ratio  was  0.15.  - 

The  disk  shapes  and  centrifugal  stress  distributions  for  the 
first-,  second-,  and  third-stage  disks  of  the  selected  turbine 
configuration  are  presented  in  Figures  11,  12,  and  13,  re¬ 
spectively. 
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TABLE  III.  SUMMARY  OF  PERTINENT  TURBINE- 
ROTOR  DATA  FOR  TURBINE  DESIGN 
A-338-5 


Turbine  Stage  Number 

Pertinent  Items 

1 

2 

2 

2 

3 

3 

Turbine  speed,  rpm  x  10”3 

38 

38 

38 

38 

38 

38 

Stage  work  split 

Equal 

Equal 

Equal 

Equal 

Equal 

Equal 

Chord  ratio,  TC / CX 

0.15 

0.15 

0.15 

0.20 

0.15 

0.20 

Number  of  blades 

40 

40 

30 

40 

40 

40 

Tip  diameter,  in. 

8.116 

8.856 

8.856 

8.856 

9.828 

9.828 

Inlet  temperature,  °F 

2400 

2076 

2076 

2076 

1752 

1752 

Blade  height,  in. 

0.74 

1.11 

1.11 

1.11 

1.60 

1.60 

Maximum  blade  stress,  ksi 

9.8 

12.6 

12.6 

12.5 

39.0 

40.0 

Disk  weight,  lb 

7.59 

7.89 

9.82 

10.70 

11.01 

16.24 

Maximum  disk  stress,  ksi 

77.0 

82.0 

77.0 

78.0 

101.0 

103.0 

Average  tangential  disk 
stress,  ksi 

44.0 

48.5 

44.0 

45.5 

61.5 

64.0 

Tip  speed,  ft  per  3ee 

1350 

1470 

1470 

1470 

1630 

1630 

Notes : 
I. 


Blade  material  for  the  first-  and  second-stage 
rotors  was  considered  to  be  beryllide  and  for 
the  third-stage  rotor  was  considered  to  be 
INCO  713  C,  as  cast. 


2.  Disk  material  for  all  three  rotor  stages  was 
considered  to  be  INCO  713 C,  as  cast. 


STRESS,  RSI 


NOTES:  1.  5PEED=38,000  RPM,  TIP  SPEED=1350  FPS 
2.  40  SLADES 
3!  TC/CX=0.15 

4.  BERYLIIDE  BLADES,  DENSITY  0.118  LB  PER  C U  IN. 

5.  DISK  MATERIAL«INCO  71 3C.  AS  CAST 

6.  INLET  TEMPERATURE-2860  °R 

7.  WEIGHT=7.59  LB 

8.  AVERAGE  TANGENTIAL  STRESS=44,000  PSI 
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0  12  3 

MUSS.  HI. 

Figure  11.  Preliminary  Disk  Shape  and 

Calculated  centrifugal  Stresses 
of  First-Stage  Turbins. 
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NOTES:  1.  SPEED=38.000  RPM.  TIP  SPEED=1470  FPS 

2.  40  BLADES 

3.  TC/CX=0.15 

4.  BERYLUDE  BLADES.  DENSITY  0.118  LB  PER  CU  IN. 

5.  DISK  MATERIAl=INCO  713C,  AS  CAST 

6.  INLET  TEMPERATURE=253fc  °R 

7.  WEIGHT=7.89  LB 

8.  AVERAGE  TANGENTIAL  STRE$S=48,500  PSI 


0  1.0  2.0  3.0 

HADISS.  IN. 


Figure  12.  Preliminary  Disk  Shape  and 

Calculated  Centrifugal  Stresses 
of  Second- Stage  Turbine. 
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STRESS,  XSI 
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NOTES:  1.  SPEED=38,000  RPM,  TIP  SPEED-1630  FPS 

2.  40  BLADES 

3.  TC/CX=0.20 

4.  DISK-BLADE  MATERIAL=INCO  713C,AS  CAST 

5.  INLET  TEMPERATURE=2212  °R 

6.  WEIGHT=16.24  LB 

7.  AVERAGE  TANGENTIAL  STRESS=64,C00  PSI 


Figure  13.  Preliminary  Disk  Shape  and 

Calculated  Centrifugal  Stresses 
of  Third-Stage  Turbine. 
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?. . 8  PRELIMINARY  THERMAL  ANALYSES  -  ROTOR  DISKS 

The  preliminary  thermal  analyses  of  the  rotor  disks  were  con¬ 
ducted  in  twc  phases.  Phase  I  studies  were  made  for  the 
first-  and  second  -stage  turbine  disks,  and  Phase  II  studies 
were  made  for  the  third-stage  turbine  disk. 

2.6^1  Phase  I  Studies 


The  first-  and  second-stage  turbine  disk  thermal  analyses 
involved  setting  up  a  sector  of  the  rotor  disks  on  an  existing 
steady-state- temperature  computer  program.  The  included  angle 
of  the  rotor  disk  sectors  was  simply  360°  divided  by  40,  the 
number  of  rotor  blades. 

Included  in  the  analyses  were  the  effect  of  thermal  radiation 
between  the  disk  and  blade  metal  surfaces  not  in  direct  con¬ 
tact,  the  surface  contact  resistance  between  the  disk  and  blade 
metal  surfaces  that  were  in  direct  contact,  and  the  changes  of 
the  metal  physical  properties  with  temperature.  Also  included 
were  the  thermal  resistance  between  the  hot  gas  and  the  blade, 
the  thermal  resistance  to  heat  flowing  from  the  disk  to  the 
cooling  air,  and  the  thermal  capacity  of  the  cooling  air. 

The  cooling  air  used  in  the  analyses  was  assumed  to  be  part 
of  the  discharge  of  a  compressor  group  with  an  overall  pres¬ 
sure  ratio  of  10:1.  It  was  assumed  that  approximately  2  per¬ 
cent  of  the  compressor  threugh-flow  would  be  required  to 
adequately  cool  the  rotor  disks.  It  was  further  assumed  that 
a  nominal  split  of  the  cooling  air  would  be  50  percent  for 
the  first  stage,  30  percent  for  the  second  stage,  and  2G  per¬ 
cent  for  the  third  stnge.  The  cooling-air  flow  path  was  the 
same  for  each  disk  and  was  assumed  to  enter  at  the  inside 
diameter  (ID)  of  the  disk,  to  flow  out  radially  with  equal 
distributions  over  both  sides  of  the  disk,  and  to  enter  the 
hot  gas  stream  at  the  outside  diameter  (OD)  of  the  disk. 

For  these  studies*-  the  blade  material  was  considered  to  be 
beryllide  and  the  disk  material  was  considered  to  be  a  nickel- 
base  alloy. 

A  number  of  computer  analyses  were  conducted  to  find  the  ef¬ 
fect  on  the  steady-state  disk  temperature  distributions  that 
variation  of  each  of  the  following  parameters  would  produce. 

1.  Cooling-air  inlet  temperature 

2.  Coolinc-air  flow  rate 

3.  Thermal  radiation  resistances 
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4. 


Hot- gas  total  temperature 


5.  Thermal  resistance  between  the  hot  gas  and  the 
blade  surface 

Hie  procedure  used  was  to  maintain  all  the  parameters  at  their 
calculated  values  (shown  in  Table  IV)  except  one,  which  was 
varied  to  determine  its  effect  on  the  pertinent  disk  tempera¬ 
tures. 

For  the  analyses,  thermal  radiation  resistances  (Item  No.  3 
above)  were  calculated  from  an  expression  involving  the  sur¬ 
face  emissivity  and  the  view  factors  between  the  surfaces. 
However,  since  the  surface  emissivity  may  vary  over  a  wide 
range  during  the  life  of  the  turbine  as  a  result  of  deposits 
on  the  metal  surfaces,  it  was  considered  important  that  the 
thermal  radiation  resistance  be  one  of  the  variable  param¬ 
eters  in  the  Phase  I  studies. 


TABLE  IV.  CALCULATED  VALUES  OF  THERMAL 

PARAMETERS  FOR  "FIRST-  AND 

SECOND-STAGE  ANALYSES 

Parameter 

First 

Stage 

Second 

Stage 

Cooling  inlet  temperature,  °r 

1160 

1160 

Cooling  flow  rate,  percent 
(total  through-flow  of 
engine,  5  lb/ sec) 

1.0 

0.6 

Radiation  resistance  factor 

| 

1.0 

1.0 

Hot-gas  total  temperature 
relative  to  the  turbine,  °R 

2620 

2290 

Effective  heat-transfer 
coefficient  between  the 
turbine  and  the  working 
fluid,  Btu/hr-fta-°F 

1760 

1162 

In  addition  to  the  analyses  that  were  made  varying  the  subject 
parameters,  analyses  were  conducted  with  all  the  parameters 
maintained  at  their  calculated  values  to  determine  the  effect 
produced  by  placing  nickel  shim  stock  between  the  mating 
surfaces  of  the  blade  and  disk. 

It  was  determined  that  with  the  shim  in  place,  the  operating 
temperature  of  the  disk  OD  is  approximately  30°  to  40°  cooler 
on  both  the  first-  and  second-stage  disks.  In  the  operational 
temperature  range  anticipated  for  these  two  turbine  disks,  a 
35°  drop  in  metal  temperature  would  allow  an  increase  in 
stress  level  of  approximately  5000  psi  with  no  change  in 
stress-rupture  life.  Also,  with  the  shim  stock  in  place,  the 
temperatures  of  the  first-  and  second-stage  blade  shanks  in¬ 
creased  approximately  30°  and  20*,  respectively,  and  the  tem¬ 
perature  difference  between  the  blade  shanks  and  the  disks 
increased  approximately  55°.  These  temperature  changes  all 
produce  desirable  effects  for  the  subject  design.  For  these 
reasons,  plus  the  added  benefit  of  having  a  relatively  soft 
material  between  the  brittle  blade  material  and  the  disk  con¬ 
tact  surfaces,  it  was  recommended  that  the  shim  be  included 
in  the  blade -attachment  final  design  and  analyses. 

Among  the  parameters  cf  Table  IV  that  were  varied  for  the 
thermal  analyses,  the  cooling-air  flow  rate  has  the  greatest 
effect  on  engine  performance  and  disk  cooling,  and  can  be 
easily  varied  in  actuality.  Therefore,  the  effect  of  this 
parameter  on  pertinent  disk  temperatures  illustrates  the 
results  of  the  Phase  I  studies.  Figures  14  and  15  show  the 
temperatures  at  various  locations  for  the  first-  and  second- 
stage  rotors,  respectively,  with  shims  assumed  in  place  ana 
as  a  function  of  the  cooling-air  flow  rate. 

2.8.2  Phase  II  Studies 

The  analytical  procedure  followed  during  the  third-stage 
thermal  analysis  was  the  same  as  that  used  for  the  Phase  I 
studies  reported  in  Paragraph  2.8.1.  The  parameters  that 
ware  varied  for  the  third-stage  analyses  included  the  fellow- 
ing: 

♦ 

% 

1.  Cooling-air  inlet  temperature 

2.  Cooling-air  flow  rate 

3.  Cooling-air  heat-transfer  coefficients 

4.  Rot-gas  total  temperature 
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Figure  14.-  Preliminary  Thermal  Analysis 
of  First-Stage  Rotor. 
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•  .  Thermal  resistance  between  tne  hot  jas  an 
the  blade  surface 

The  procedure  used  was  to  maintain  all  the  parameters  at  their 
calculated  values  (shown  in  Table  V)  except  one,  which  was 
varied  to  determine  its  effect  on  the  pertinent  disk  tempera¬ 
tures. 


TABLE  V.  CALCULATED  VALUES  OF 

THERMAL  PARAMETERS  FOR 
THIRD- STAGE  ANALYSES 


Cooling  inlet  temperature,  °R 

Cooling  flow  rate,  percent  (total 
through-flow  of  engine,  6.0  lb/ 
sec) 

Factor  applied  to  the  calculated 
heat-transfer  coefficients 
between  the  cooling  gas  and 
disk  sides 

Hot  gas  total  temperature  relative 
to  the  turbine,  °R 

Effective  heat-transfer  coefficient 
between  the  turbine  and  the 
working  fluid,  Btu/hr-ft2-°F 


1160 

0.4 

1.0 

1970 

676 


The  third-stage  turbine  disk  and  blades  were  designed  as  an 
integrally  cast  one-piece  rotor  of  IN-100  material.  The 
shape  of  the  disk  was  patterned  after  tha  TPE331  turbine 
rotors,  with  the  rim  slotted  between  every  blade  to  reduce 
the  calculated  tangential  stresses  due  to  thermal  effects 
near  the  OD  of  the  hub.  A  hole  is  drilled  through  the  rotor 
at  the  bottom  of  each  slot  to  reduce  the  stress  concentration 
that  would  normally  exist  at  the  bottom  of  the  slot. 

Again,  as  in  the  Phase  I  studies,  the  effect  of  the  cooling- 
air  flow-rate  parameter  illustrates  the  results  of  the  thermal 
study.  Figure  16  shows  the  temperature  at  various  locations 
for  the  third-stage  rotor  as  a  function  of  the  cooling-air 
flow  rate. 
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Figure  16.  Preliminary  Thermal  Analysis 
of  Third-Stage  Rotor. 


40 


2.9  PRELIMINARY  C ENTRI i't’GAL-  AND  THERMAL- STRESS  ANALYSIS 


An  clastic  stress  analysis  was  conducted  to  calculate  the 
combined  thermal-  and  centrifugal-stress  distribution  for 
each  turbine  disk.  The  studies  revealed  a  major  stress  prob¬ 
lem  in  each  of  the  disks,  in  that  the  combined  thermal  and 
centrifugal  stresses  that  exist  near  the  disk  bore  of  each 
stage  (tangential  stress — tensile)  and  near  the  blade  base 
of  the  first-stage  disk  (tangential  stress — compressive)  are 
beyond  the  yield  strength  of  the  material.  Figures  17,  18, 
and  19  present  the  results  of  this  analysis.  Shown  on  these 
figures  are  the  combined  thermal-  and  centrifugal-stress 
distributions,  the  calculated  temperature  distributions,  and 
the  stress  distributions  due  to  centrifugal  effects  alone. 

The  thermal  stresses  represented  a  problem  that  required 
solution.  The  considerable  magnitude  of  the  thermal 
stresses  was  the  result  of  the  high  temperature  and  the 
severity  of  the  radial  temperature  gradients  of  the  rotor 
disks.  The  problem,  therefore,  was  to  define  a  cooling-air 
flow  path,  or  scheme,  that  would  provide  a  more  effective 
means  of  disk  cooling  and  produce  a  more  favorable  tempera¬ 
ture  gradient  in  the  disk.  As  described  previously,  the 
cooling  air  was  assumed  to  enter  at  the  ID  of  the  disk,  to 
flow  out  radially  with  equal  distributions  over  both  sides 
of  the  disk,  and  to  enter  the  hot-gas  stream  at  the  OD  of 
the  disk. 


Thermal  studies  on  the  first-stage  disk  were  initiated  to 
determine  a  cooling  scheme  that  would  reduce  the  disk  tem¬ 
perature  levels  and  thermal  gradients  without  increasing  the 
flow  rate  of  cooling  air.  As  reported  in  Paragraph  3. 3. 1.3, 
these  efforts  were  successful  and  the  resulting  disk 
thermal  stresses  were  within  design  life  limits. 


With  the  exception  of  the  centrifugal  stresses  near  the  bore 
of  the  third-stage  disk  (tangential  stress),  the  centrifugal 
stresses  in  tho  disks  are  well  below  th'">  yield  ctror'gi-h 
of  the  material. 

With  completion  of  the  centrifugal-  and  thermal-stress  analy¬ 
ses,  the  preliminary  design  studies  were  concluded.  Efforts 
were  then  initiated  to  complete  the  final  aerody  namic  aesign 
and  to  conduct  thermal  and  stress  analyses  on  the  anticipated 
mechanical  configurations  of  the  first-stage  turbine  com¬ 
ponents. 
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NOTES: 

1.  SPEED  =  38,000  RPM,  TIP  SPEED=1350  FPS 

2.  40  BLADES,  TC/CX=0.15 

3.  TURBINE  DISK  MATERIAL-INCO  713C,AS  CAST 

4.  INLET  TEMPER  ATURE=2860  °F 

5.  WEIGHT=7.59  LB 

6.  AVERAGE  TANGENTIAL  STRESS=44,000  PSI 


0  1.0  2.0  3.0 

RADIUS,  IN. 

- COMBINED  THERMAL  AND  CENTRIFUGAL  STRESS 

- CENTRIFUGAL  STRESS  ONLY 


Figure  17.  Preliminary  Centrifugal-  and 

Thermal-Stress  Distributions  of 
First-Stage  Rotor  Disk. 
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STRESS,  KSI 


NOTES 

1.  SPEED  =  3S,000  RPM.  TIP  SPEED=I630  FPS 

2.  40  8LADES,  TC/CX=0.20 

3.  TURBINE  WHEEL  MATERIALUNCO  713C.AS  CAST 

4.  INLET  TEMPERATURES 21 2  °R 

5.  WEIGHT=16.24  LB 

6.  AVERAGE  TANGENTIAL  STRESS=64,000  PSI 


RADIUS,  IN. 

COMBINED  THERMAL  AND  CENTRIFUGAL  STRESS 
CENTRIFUGAL  STRESS  ONLY 


Figure  19.  Preliminary  Centrifugal-  and 

Thermal-Stress  Distributions  of 
Third-Stage  Rotor  Disk. 
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3.1  GENERAL 

Upon  completion  of  the  preliminary  design  studies,  calcula¬ 
tions  were  made  to  determine  the  preliminary  margins  of 
safety  for  Turbine  Design  A-338-5.  These  calculations  were 
based  on  material-property  estimates  and  data  acquired  from 
preliminary  material-properties  testing  and  on  stress  analy¬ 
ses  of  the  preliminary  design  studies.  Table  VI  presents  the 
calculated  stress  levels,  estimated  material  strengths,  and 
the  calculated  margins  of  safety. 

The  final  design  activities  were  initiated  with  aerodynamic 
analyses  to  establish  the  first-stage  turbine  stator-vane 
and  rotor-blade  detailed  shapes  and  to  refine  the  overall 
gas-flow  path  of  Turbine  Design  A-338-5.  Mechanical  studies 
were  then  conducted  to  determine  the  most  advantageous 
structural  configuration  of  the  first-stage  turbine  com¬ 
ponents  without  affecting  the  established  airfoil  shapes, 
and  to  analyze  the  final  component  configurations  for  two 
potential  failure  modes — stress  rupture  and  low-cycle  fa¬ 
tigue.  Results  of  these  activities  are  reported  on  the 
following  pages. 

3.2  AERODYNAMIC  FINAL  DESIGN  ACTIVITIES 

Aerodynamic  design  studies  were  conducted  to  refine  the 
overall  gas-flow  path  of  Turbine  Design  A-338-5  and  to  define 
the  detailed  stater  vane  and  rotor  blade  shapes  for  the  first- 
stage  turbine  components.  Further  aerodynamic  studies  were 
to  be  performed  when  the  mechanical  studies  were  completed 
and  anticipated  results  of  the  cascade  testing  (reported  in 
Volume  I)  were  known. 

3.2.1  Gas-Flow- Path  Refinement 


New  vector  diagrams  were  made  for  the  A-338-5  three-stage 
t  nine  design  to  reflect  the  lower  loss  coefficients  ob- 
ta_ned  from  a  blading  efficiency  analysis.  Representative 
vector  diagrams  for  the  first-stage  turbine  stator  and  rotor 
are  shown  in  Figure  20. 

The  hub  and  tip  diameters  of  the  stators  and  rotors  of  each, 
of  the  three  stages  were  revised  slightly  so  that  for  normal 
machining  and  assembly  tolerances,  each  succeeding  stator  an  1 
rotor  would  be  smaller  in  diameter .  This  ensured  that  a 
"step-up"  ir.  flow  would  not  occur  through  the  stages  that 
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TABLE  VI.  PRELIMINARY  CALCULATIONS  FOR 

MARGINS  OF  SAFETY 

Location 

Stress 

(psi) 

Temn 

(°F) 

Strength 

(psi) 

Margin 

of 

Safety 

(M.S.) 

First-Stage  Turbine 

Blade  root  (beryl lice) 

10,000 

2140 

23,000  SR 

1.3 

Blade  attachment  (beryl lide) 

17,500 

1900 

37,000  SR 

1.1 

Disk  attachment  (INCO  713C) 

17,500 

1450 

50,000  SR 

1.85 

Disk  average  tangential 
stress  (INCO  713C) 

44,000 

1200 

85,000  YS 

0.93 

Second-Stage  Turbine 

31ade  root  (beryllide) 

13,000 

1820 

43,000  SR 

2.31 

Blade  attachment  (beryllide) 

28,000 

1640 

50,000  SR 

0.78 

Disk  attachment  (INCO  713C) 

28,000 

1550 

70,000  SR 

1.5 

Disk  average  tangential 
(INCO  713C) 

48,500 

1200 

85,000  YS 

0.75 

Third-Stage  Turbine 

Blade  root  (INCO  713C) 

3  9,000 

1500 

47,000  SR 

0.20 

Disk  average  tangential 
(INCO  713C) 

64,000 

1200 

85,000  YS 

0.32 

NOTES : 

2.  YS  =  material  yield 

strength  at  temperature. 

3.  SR  =  material  stress 
for  200-hour  11 

-rupture 

fe. 

strength  at  temperature 

4.  The  strength  levels  chosen  for  beryllide  parts  are  based 
on  the  anticipated  strength  levels  that  would  be  achieved 
in  the  final  beryllide  alloy.  There  is  sufficient  range 
in  the  present  design  to  allow  for  as  much  as  a  50-per¬ 
cent  reduction  in  stress  level  in  the  beryllide  parts  anc 
still  maintain  a  satisfactory  margin  of  safety  in  the 
disk  portions  of  the  turbines. 

5.  Rotational  speed  =  38, 

000  rpm. 
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could  affect  turbine  efficiency.  Table  VII  reflects  the  old 
and  the  new  revised  turbine  stage  dimensions  of  Turbine 
Design  A-338-5.  The  meridional  view  of  the  first  stage  is 
shown  in  Figure  21,  with  the  radii  indicated  for  the  local 
vector  diagrams  and  detailed  vane  and  blade  shape  studies. 


3.2.2  Detailed  Stator-Vane  Shape 

The  final  detailed  stator-vane  shape  is  shown  in  Figures  22 
and  23  for  radii  of  3.5694  and  3.9694  inches.  These  radii 
correspond  to  approximately  one-quarter  and  three-quarters 
of  the  vane  height,  respectively.  The  stacking  line  for 
these  vane  sections  was  chosen  at  the  center  of  the  throat 
and  was  used  as  the  origin  of  coordinates.  Coordinates  for 
vane  sections  at  other  radii  were  obtained  from  a  linear 
interpolation  or  extrapolation  by  a  computer  program. 

Also  shown  on  Figures  22  and  23  are  the  corresponding  surface 
velocity  distributions  in  terms  of  velocity  ratio,  v/a  « ,  as 

a  function  of  the  percentage  of  the  axial  chord.  As  antici¬ 
pated,  the  surface  velocity  distributions  show  high  reaction 
(large  acceleration  of  the  fluid)  and  minimum  diffusion  of 
the  surface  velocity,  and  therefore,  represent  an  acceptable 
aerodynamic  stator-vane  design. 


3.2.3  Detailed  Rotor  Blade  Shape 


The  final  detailed  rotor  blade  shape  was  based  on  analyses 
conducted  at  five  cross  sections.  These  five  cross  sections 
were  taken  at  the  following  radii  and  corresponding  blade 
positions : 


1. 

3.4204 

inches 

2. 

3.5694 

inches 

3. 

3.7694 

inches 

4. 

3 . 9694 

inches 

5. 

4.1057 

inches 

blade  root 

one-quarter  of  blade  height 
mean  blade  height 
three-quarters  of  blade  height 
blade  tip 


The  stacking  line  for  these  blade  sections  is  the  center  of 
gravity  of  each  section.  Coordinates  for  blade  sections  at 
other  radii  were  obtained  from  a  linear  interpolation  or 
extrapolation  by  a  computer  program. 
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TABLE 

VII.  TURBINE  A- 3 38- 5 

STAGE  DIMENSIONS 

Old  Dimensions 

New  Dimensions 

Hub  Radius 

Tip  Radius 

Hub  Radius 

Tip  Radius 

Stages 

(in.) 

(in.) 

(in.) 

(in. ) 

First 

Stator  • 

Inlet 

3.318 

4.058 

3.4204 

4.1057 

Exit 

3.318 

- 

3.4204 

4.1057 

Rotor 

Inlet 

3.318 

4.058 

3.4102 

4.1452 

Exit 

3.318 

- 

3.4102 

4.1452 

Second 

Stator 

Inlet 

3.318 

4.058 

3.4000 

4.1452 

Exit, 

3.318 

4.428 

3.4000 

4.4405 

Rotor 

Inlet 

3.318 

4.428 

3.3898 

4.4803 

Exit 

3.318 

- 

3.3898 

4.4303 

Third 

Stator 

Inlet 

3.318 

4.428 

3.3796 

4.4803 

Exit 

3.318 

4.914 

3.3796 

4.9085 

Rotor 

Inlet 

Exit 


3.314 

3.314 


4.914 


3.3694 

3.3694 


4.9496 

4.9496 
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Figure  21.  Meridional  View  of  First;- 5'.. age  Turbine 


VELOCITY  RATIO  ,  V  'a'  PERCENT  OF  VANE 

AXAIL  CHORD  »C 


The  blade  shape  is  shown  in  Figures  24  and  25  for  radii  of 
3.5694  and  4.1057  inches.  Also  shown  in  these  figures  are 
the  corresponding  surface  velocity  distributions  in  terms 
of  velocity  ratio,  W/a^,  as  a  function  of  the  percentage 
of  the  axial  chord.  Evaluation  of  these  data  shows  that 
only  slight  suction  surface  diffusion  is  present  and  the 
overall  blade  reaction  increases  rapidly  with  an  increasing 
radius  and,  therefore,  the  rotor  blade  shape  represents  an 
acceptable  aerodynamic  rotor  blade  design. 

3.3  MECHANICAL  STUDIES 

With  the  detailed  aerodynamic  shapes  established  for  the 
first-stage  turbine  stator  vane  and  rotor  blade,  mechanical 
studies  were  initiated  to  determine  the  most  advantageous 
structural  configuration  of  the  first-stage  turbine  com¬ 
ponents  without  affecting  the  established  airfoil  shapes. 
Included  in  the  mechanical  studies  were  detailed  thermal  and 
thermal-stress  analyses,  stator-vane  and  rotor-blade  attach¬ 
ment  studies,  and  photoelastic  and  thermal-conduction  test 
evaluations. 

The  mechanical  design  goal,  as  stated  in  the  preliminary 
design  studies,  was  a  200-hour  stress-rupture  life  at  2400°F 
TIT.  The  final  design  mechanical  studies  included  sufficient 
stress  analyses  to  facilitate  evaluation  of  two  potential 
failure  modes — stress  rupture  and  low-cycle  fatigue.  For 
the  stress-rupture  evaluation,  it  was  planned  that  the  per¬ 
tinent  calculated  stresses  from  the  mechanical  studies  would 
be  compared  to  the  stress-rupture  properties  of  the  turbine 
material,  which  were  to  be  determined  in  the  concurrent 
turbine-materials  research  program  reported  in  Volume  I. 

For  the  low-cycle  fatigue  evaluation,  it  was  planned  that 
the  pertinent  calculated  stresses  from,  the  mechanical  studies 
would  be  compared  to  those  determined  by  similar  analyses  in 
operational  gas-turbine  engines  of  comparable  size  and  design. 
No  attempt  was  planned  at  this  time  to  obtain  low-cycle 
fatigue  properties  for  the  beryllide  material  because  of 
the  complexities  involved  in  this  type  of  testing.  However, 
this  task  was  anticipated  to  be  accomplished  later  in  the 
program,  depending  upon  the  degree  of  success  that  was  real¬ 
ized  with  the  beryllide  material. 

At  this  time  in  the  program,  the  material  properties  of 
IN-100  (rotor-disk  material)  had  been  established,  mate¬ 
rial  property  testing  of  Composition  1  (stator-vane 
material)  had  been  initiated,  the  stator- vane  cascade  test 
rig  was  being  developed,  and  Composition  2  (rotor-blade 
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PFRCENT  OF  BLADE  AXIAL  CHORD  C 

x 

First-Stage  Rotor  Blade  Shape 
ana  Corresponding  Surface 
Velocity  Distribution  for 
&  —  3.5694  Inches. 


Figure  24. 
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First-Stage  Rotor  Blade 
Shape  and  Corresponding 
Surfpce  Velocity  Distribution 
for  R  =  4.1057  Inches. 


Figure  25 


! 


material)  material-development  studies  had  not  yet  been  ini¬ 
tiated  (reported  in  Volume  I).  Therefore,  it  was  necessary 
to  estimate  the  physical  properties  of  the  beryllide  material . 
The  pertinent  physical  properties  that  were  estimated  for 
use  in  these  mechanical  studies  are  shown  in  Table  VIII.  For 
information  purposes.  Table  VIII  also  includes  measured 
values  of  the  physical  properties  that  were  obtained  at  a 
later  date  from  the  materials  research  program. 


TABLE  VIII.  COMPARISON  OF  BERYLLIDE  PHYSICAL 
PROPERTIES  ESTIMATED  AND 
MEASURED  AT  800°F 


Physical  Property 

Estimated 

Values 

Measured 

Values 

Modulus  of  elasticity, 
psi 

43.9  x  10s 

42.3  x  10s 

Coefficient  of  thermal 
expansion, 
in. /in.  -  °F 

8.83  x  1G~3 

8.20  x  10"6 

Density,  lb/cu  in. 

0.  114 

0.103 

Specific  heat, 

Btu/lb-  CF 

0.340 

0.312 

Thermal  conductivity, 

Bt  u/hr- f t 2 - °  F/ in . 

222 

165 

3.3.1  Rotor — Mechanical  Design  Studies 

The  mechanical  design  studies  that  were  conducted  to  deter¬ 
mine  and  evaluate  structural  configurations  of  the  rotor 
blade,  rotor  disk/blade  attachment,  and  rotor  disk  of  the 
first- stage  turbine  are  reported  ir.  the  following  sub- 
paragraphs. 
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3 . 3 . 1  - 1 


Rotor  Blade 


The  first  mechanical  configuration  of  the  aerodynamic  blade- 
design  to  be  evaluated  was  that  of  a  solid  blade.  Initial 
studies  of  the  solid  tlade  revealed  that  the  centrifugal 
blade  load  imposed  on  the  disk  was  approximately  2680  pounds 
per  blade.  This  value  v/as  nearly  twice  the  estimated  blade 
loading  used  for  the  preliminary-design  disk  shape  studies 
and  was  considered  to  be  prohibitive  to  the  anticipated  rotor 
disk/biade  attachment.  To  reduce  the  magnitude  of  the  blade 
loading,  studies  were  initiated  to  design  a  radial  hole  in 
the  solid  blade  and,  thus,  reduce  the  blade  mass.  Analyses 
were  conducted  to  determine  the  most  desirable  hole  size  and 
location  to  minimize  the  thermal  and  centrifugal  loading 
stresses. 

In  determining  the  contour  of  the  radial  hole,  a  minimum 
blade  wall  thickness  of  0.025  inch  was  selected.  Tnis  dimen¬ 
sion  was  located  near  the  minimum  value  of  heat-transfer 
coefficient  on  the  pressure  side  of  the  blade,  since  «:he 
wall  thicknesses  around  the  radial  hole  were  determined  with 
the  objective  of  minimizing  the  thermal  gradients  occurring 
over  the  cross  section  at  any  blade  section.  Because  of  the 
blade  taper,  from  tip  to  root,  the  minimum  wall  thickness 
varied  from  0.025  inch  at  the  tip  end  tc  0.032  inch  at  the 
root  end  of  the  hole.  The  radial  hole  reduced  the  maos  of  the 
airfoil  portion  of  the  blade  by  41  percent.  The  solid  blade 
centrifugal  load  of  2680  pounds  per  blade  was  reduced  to  1640 
pounds  per  blade  after  adding  the  radial  hole.  The  shape  of 
the  contoured  radial  hole  in  the  blade  is  shown  in  Figure  26. 
Also  shown  in  this . figure  is  the  calculated  static-pressure 
distribution  around  the  blade  at  the  cross  section  located 
at  a  radius  of  3.7694  (mean  blade  height) .  From  the  calcu¬ 
lated  pressure  distributions,  the  air  loads  'were  determined 
for  the  disk/blade  attachment  studies. 

Detailed  thermal  and  thermal-stress  analyses  were  conducted 
on  the  rotor  blade.  For  comparison  purposes,  the  studios 
included  both  hollow  and  solid  blades.  The  analyses  involved 
determination  of  the  steady-state  heat-transfer-coef f iciont 
distribution  and  relative  gas-temperature  distribution  around 
the  blade  (these  data  were  available  for  the  radial  hole 
contour  studies),  the  transient  temperature  distributions, 
and  the  corresponding  thermal-stress  distributions  of  the 
blade.  Variations  in  gas  temperature  and  heat- transfer-  coef¬ 
ficients  with  time  were  included  in  the  analyses,  as  were  the 
blade  material- property  variations  with  temperature.  The 
time-dependent  properties  were  the  result  of  calculating  the 
TIT  variation  with  time  based  on  existing  experimental  start 
data  from  three-stage  turbines  of  lower  operating  TIT  ranges 
(see  Figure  27'*  . 
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Other  assumptions  made  for  the  thermal-stress  calculations  of 
the  blade  were: 

1.  There  were  no  thermal  gradients  in  the  radial 
direction. 

2.  All  stresses  were  due  to  thermal  gradients  in 
the  remaining  two  directions  {no  external  loads) . 

The  following  conditions  of  equilibrium  were  satisfied  at 
the  cross  section  analyzed: 


J  a  dA  =  0 


ox  dA  =  0 


J  ay  dA  =  o 

where  x  and  y  are  any  axes  in  the  plane  of  the  cross  section 
and  A  is  the  cross-sectional  area  of  the  blade  section. 

The  heat-transfer-coefficient  distribution  and  the  relative 
gas- temperature  distribution  around  the  blade  are  shown  in 
Figure  28  at  steady  state.  A  large  variation  in  heat- 
transfer  coefficient  exists  across  the  thick  portion  of  the 
blade,  which  results  in  large  thermal  gradients  in  the  blade 
during  start-up.  Also,  the  large  variation  in  blade  thick- 
ness  from  the  middle  of  the  blade  cross  section  to  the 
trailing  edge  increases  the  magnitudes  of  the  thermal  grad¬ 
ients  that  occur  and  the  corresponding  thermal  stress  levels. 

Figures  29  and  30  show  the  results  of  the  transient  tempera¬ 
ture  and  thermal-stress  analyses  of  the  rotor  blade  cross 
section  at  a  radius  of  4.1057  inches  (blade  tip).  The  blade 
tip  cross  section  was  selected  for  these  analyses,  since 
it  was  the  farthest  removed  from  the  blade  root  and,  there¬ 
fore,  the  least  likely  to  be  affected  by  radial  conduction 
of  heat  into  the  disk/blade  attachment.  Experience  has 
shown  that  radial  conduction  of  heat  from  the  blade  cross 
sections  near  the  root  into  the  disk/blade  attachment  region 
tends  to  produce  thermal  gradients  that  are  less  severe  than 
those  predicted  by  a  two-dimensional  heat-transfer  analysis 
that  neglects  radial  conduction.  As  shown  in  Figure  29, 
making  the  blade  hollow  caused  an  approximate  25-percent 
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Figure  28.  Steady-State  Heat-Transfer-Coefficient 

Distribution  and  Relative  Gas— Temperature 
Distribution  Around  the  Blade. 


NOTES: 


1  SEE  FIGURE  27  FOR  FIT 
VARIATION  VERSUS  TIME 

3.  STRESS  AND  TEMRERATURE 
ROINTS  RIOTTED  ARE  IOCATED 
NEAR  ROINTS  OR  MAXIMUM 
STRESS 

3.  TIME  MEASURED  FROM 
COMRUSTOR  IIGHT.QFF 


Figure  29.  Transient  Thermal  Study  of 
Rotor  Blade  Section  at 
R  =  4.1057  Inches  (Blade  Tip). 


reduction  in  the  peak  positive  thermal  stress  level  and  an 
approximate  16-percent  increase  in  the  peak  compressive 
stress  level.  The  solid  blade  material  temperature  at  the 
time  and  place  these  peak  stresses  occur  is  approximately 
1100°F  and  700°F  at  stress  levels  of  minus  90,000  psi  and 
plus  75,000  psi,  respectively.  For  the  hollow  blade,  the 
temperatures  are  approximately  the  same,  and  the  stress 
levels  are  minus  105,000  psi  and  plus  55,000  psi,  respec¬ 
tively. 

These  transient  thermal  stresses  are  of  a  magnitude  (80,000 
to  120,000  psi)  similar  to  that  calculated  in  operational 
gas-turbine  engines  of  comparable  size  and  design.  No  further 
attempt  was  made  at  this  time  to  reduce  the  transient  thermal 
stresses  or  to  perform  a  low-cycle- fatigue  analysis  on  the 
subject  rotor  blade  configurations.  As  previously  stated, 
low-cycle- fatigue  analyses  were  to  be  conducted  in  a  later 
phase  of  the  program  provided  that  the  beryllide  material 
successfully  satisfied  the  stress-rupture  requirements  of 
this  phase  of  the  final  design. 

The  transient  thermal  stress  analyses  (Figures  29  and  30) 
were  performed  to  investigate  whether  the  estimated  physical 
properties  of  the  beryllide  turbine  material  would  result  in 
excessive  thermal  stress  levels  as  compared  to  those  calcu¬ 
lated  in  conventional  turbine  material.  Since  nc  stresses 
encountered  were  significantly  greater  than  those  occurring 
in  conventional  turbine  designs,  and  beryllide  material 
properties  were  not  yet  available  from  the  materials  test 
program,  the  hollow-blade  configuration  was  judged  to  be 
an  acceptable  mechanical  design. 

3. 3. 1.2  Rotor  Disk/Blade  Attachment 

The  rotor  disk/blade  attachment  design  objective  was  to 
determine  a  configuration  in  which  the  calculated  maximum 
stress  in  both  the  disk  and  the  blade  portions  did  not 
exceed  40,000  psi  at  operating  speed.  This  stress  level 
was  again  based  on  preliminary  materials  testing  with  the 
beryllide  compositions. 

It  was  evident  that  a  dovetail  attachment  design  (single 
lobe)  required  considerably  less  ductility  in  the  disk/blade 
material  combination  to  obtain  an  even  distribution  of  the 
loading  than  that  required  for  a  fir-tree  attachment  design 
(double  lobe) .  For  this  reason,  the  disk/blade  attachment 
studies  were  conducted  initially  for  a  dovetail  attachment 
design. 
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3 . 3 . 1 . 2 . 1  Dovetail  Attachment  Design  and  Analyses 

Figure  31  shows  the  general  configuration,  stress  point  lo¬ 
cations  and  pertinent  dimension  locations  of  the  dovetail 
attachment  design  that  was  analyzed.  In  the  design  of  dove¬ 
tail  blade  attachments,  the  critical  location  for  stress  in 
the  blade  is  normally  where  the  necked  section  joins  the 
tang  or  hook  of  the  dovetail.  This  section  is  located  be¬ 
tween  points  B  and  D  in  Figure  31.  Because  of  the  non¬ 
uniformity  of  load  distribution  and  neck  stress  distribution, 
the  fillet  stress  at  this  section  can  normally  be  expected 
to  be  highest  on  one  of  the  four  external  corners  labeled 
1,  2,  3,  and  4  in  Section  A -A  of  Figure  31.  The  corner 
having  the  critical  stress  will  depend  on  the  magnitude  and 
direction  of  the  leads.  Similar  circumstances  exist  for  the 
disk  portion  of  the  attachment. 

The  fillet  stress,  is  a  measure  of  nearness  to  failure 

and  is  actually  a  combination  of  the  neck  stress,  7  ,  and  the 
tang  stress,  -t^..  Points  B  and  D  in  Figure  31  are  the  loca¬ 
tions  of  the  blade  neck  stress  and  tang  stress,  respectively. 
The  equation  that  relates  to  and  can  be  stated  as: 


where 


’f  =  °a  T-t 


T  = 


1  t  c  \aa/3t) 


(23) 


c  = 


2 


and  \  is  the  angular  separation  (in  degrees)  between  the  tang 
stress  and  neck  stress  locations  in  Figure  31.  A  value  of 
44  for  N  was  used  in  the  analyses. 

For  the  disk  stress  analysis,  the  same  analytical  approach 
was  used.  The  maximum  stress  was  calculated  to  occur  at  one 
of  the  four  external  corners  labeled  on  Section  C-C  located 
between  Points  E  and  F  in  Figure  31.  The  fillet  stress  at 
each  corner  was  determined  from  the  same  relationship  between 
tang  and  neck  stresses  as  used  in  the  blade.  In  Figure  31, 
Point  E  is  the  location  of  the  disk  tang  stress,  and  Point  F 
is  the  location  of  the  disk  neck  stress. 
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Figure  31.  General  Configuration,  Stress  Point 
Locations,  and  Pertinent  Dimension 
Locations  of  Dovetail  Attachment. 
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The  location  of  Point  E  in  the  disk  and  Point  D  in  the  blade, 
where  the  tang  stresses  are  calculated,  depends  on  the  value 
assumed  for  the  angle,  y.  For  single-tang  dovetail  shapes, 
which  lack  the  stiffening  of  additional  material  under  them 
and  which  have  relatively  steep  flank  angles  (23),  the  point 
of  maximum  stress  may  be  very  close  to  the  tangency  point. 
Thus,  ¥  may  be  taken  as  low  as  5°.  For  these  analyses,  angle 
y  was  assumed  equal  to  10°  for  both  the  blade  and  disk. 

In  the  stresses  that  were  tabulated,  a  nominal  and  a  maximum 
fillet  stress  occurred.  The  difference  between  these  two 
is  that  the  bending  and  tensile  stresses  that  make  up  the 
fillet  stress  have  been  combined  with  the  appropriate  stress- 
concentration  factors  to  give  the  maximum  calculated  fillet 
stress  at  the  locations  shewn.  This  information  was  con¬ 
sidered  especially  pertinent  for  the  blade  material,  since 
the  ductility  and  the  susceptibility  to  brittle  failure  of 
the  blade  material  were  uncertain  at  this  lime. 

The  stress  analyses  conducted  for  the  dovetail  attachment 
consisted  of  studies  of  eight  different  configurations  with 
use  of  an  existing  computer  program.  Included  in  the  program 
were  axial  and  tangential  shear  loads  acting  on  the  blade, 
as  well  as  moments  resulting  from  axial  and  tangential  gas 
loads  and  bending  loads  in  the  blade.  Table  IX  presents 
a  summary  of  the  pertinent  dimensions  (defined  and  listed 
below)  or  Configuration  No.  1  and  also  shows  the  results  of 
the  stress  analysis.  The  pertinent  dimensions  used  for 
Configuration  No.  1  were  established  as  the  standard  values 
and  were  used  as  base-line  references  for  the  other  con¬ 
figurations. 

Analyses  for  each  succeeding  configuration  were  conducted 
with  one  of  the  pertinent  dimensions  varied  while  the  re¬ 
maining  pertinent,  dimensions  were  maintained  constant  at  the 
established  base-line  values.  In  this  manner,  the  effect 
that  each  of  the  pertinent  dimensions  had  on  the  calculated 
stress  values  of  the  attachment  was  evaluated.  The  pertinent 
dimensions  whose  values  were  varied  in  the  stress  studies 
were: 

1.  The  flank  angle,  29 

2.  The  blade  neck  thickness,  a 

3.  The  disk  neck  thickness,  d 

4.  The  fillet  radii,  Rg,  Rg 

5.  The  flank  lengths,  La, 
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TABLE  IX.  PERTINENT  DIMENSIONS  AND  CALCULATED 

STRESSES  OF  DOVETAIL  DISK/ BLADE  ATTACHMENT, 
_ CONFIGURATION  NO.  1  . __ 

(Refer  to  Figure  31  for  Locations) 
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With  completion  of  the  stress  analysis  on  the  seventh  con¬ 
figuration,  a  study  was  made  to  combine  the  optimum  values 
of  the  pertinent  dimensions  determined  by  the  seven  studies 
into  one  design.  Thus,  this  design  represented  the  optimum 
configuration  of  the  dovetail  attachment.  Table  X  presents 
the  dimensions  and  stresses  of  the  optimum  configuration. 

The  stress  objective  (40,000  psi)  for  the  design  was  not 
realized.  At  this  time,  efforts  were  directed  toward  design 
and  stress  analyses  of  a  fir-tree  disk-blade  attachment. 

3 . 3 . 1 . 2 . 2  Fir-Tree  Attachment  Design  and  Analyses 

Figure  32  shows  the  general  configuration,  stress-point  lo¬ 
cations,  and  pertinent  dimension  locations  of  the  fir-tree 
attachment  design  that  was  analyzed.  As  discussed  in  the 
dovetail  attachment  analyses,  the  critical  location  for 
stress  in  the  dovetail  design  is  in  the  fillet  area  of  the 
necked  sections  between  the  lobes.  This  is  also  true  for 
the  fir-tree  design.  The  fillet  stress,  o>,  of  the  fir- 

tree  design  was  calculated  in  the  same  manner  as  that  in 
the  dovetail  analyses. 

Twenty-five  fir-tree  configurations  were  analyzed  with  use 
of  a  computer  program.  Included  in  the  program  were  axial 
and  tangential  shear  loads  acting  on  the  blade,  as  well  as 
moments  resulting  from  axial  and  tangential  gas  loads  and 
bending  loads  in  the  blade.  Data  from  the  initial  con¬ 
figuration  (nominal  pertinent  dimensions  and  angles  for 
base-line  references)  and  from  the  optimized  configuration 
analyses  are  shown  in  Tables  XI  and  XII,  respectively.  The 
analyses  were  conducted  in  a  manner  similar  to  that  of  the 
dovetail  analyses,  where  all  pertinent  dimensions  and  angles 
were  maintained  at  the  base-line  values  except  one,  which 
was  varied.  The  effect  of  each  variable  was  evaluated,  and 
the  pertinent  dimension  and  angle  values  determined  to  be 
most  advantageous  were  selected  for  the  optimized  configura¬ 
tion.  The  pertinent  dimensions  and  angles  that  were  varied 
during  the  analyses  are  listed  in  the  referenced  tables. 

As  can  be  seen  in  Table  XII,  in  the  recommended  attachment 
design,  the  maximum  blade- fillet  stresses  vary  between 
31,500  psi  and  40,200  psi,  while  in  the  disk  the  nominal 
fillet  stresses  vary  between  30,000  psi  and  44,100  psi.  A 
maximum  disk  temperature  of  1500 °F  was  estimated  for  the 
point  where  the  stresses  exceed  40,000  psi  m  the  disk 
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TABLE  X.  PERTINENT  DIMENSIONS  AND  CALCULATED  STRESSES 

OF  THE  OPTIMIZED  DOVETAIL  DISK/BLADE  ATTACHMENT 
_ DESIGN,  CONFIGURATION  NO.  8 
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(NOT  TO  SCALE) 


Figure  32.  General  Configuration,  Stress  Point 
Locations,  and  Pertinent  Dimension 
Locations  of  Fir-Tree  Attachment. 


TABLE  XI .  PERTINENT  DI KEN. SION'S  AND  CALCUTATFl' 
STRESSED  OF  FIR-TREE  DISK.  BLADE 
ATTACHMENT,  CONFIGURATION  NO.  1 
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(based  on  the  results  of  preliminary  temperature  studies  that 
were  being  conducted  on  the  disk  to  determine  the  optimum 
impingement  cooling  scheme) .  The  minimum  200-hour  stress-to- 
rupture  value  for  IN-100  of  1500 °F  is  60,000  psi.  Thus,  a 
disk  stress  of  44,100  psi  at  this  point  was  acceptable. 

During  the  time  period  of  the  disk /blade  attachment  studies, 
a  program  was  being  prepared  to  experimentally  evaluate  the 
calculated  stress  distribution  occurring  in  the  disk/blade 
attachment  by  using  the  techniques  of  photoelasticity.  This 
program,  reported  in  Appendix  III  of  this  report,  consisted 
of  three  major-  parts: 

1.  Design  and  fabrication  of  a  photoelastic  test  rig. 

2.  Design  and  fabrication  of  the  experimental  attach¬ 
ment  model  to  be  tested. 

3.  Preparation  of  test  procedures,  conduction  of  the 
test,  and  analysis  of  the  data. 

Results  from  the  photoelastic  testing  indicated  that  the  fir- 
tree  attachment  stresses  were  somewhat  lower  than  the  calcu¬ 
lated  stresses  shown  in  Table  XII.  There  were  a  number  of 
differences  between  the  fir  tree  modeled  in  the  photoelastic 
tests  and  the  one  analyzed  in  Table  XII.  These  differences 
include: 

1.  An  inability  to  model  disk  curvature  in  the  photo¬ 
elastic  test  rig. 

2.  In  the  photcelastic  tests,  no  distinction  was  made 
between  the  blade  centrifugal  loading  by  itself  and 
the  combined  centrifugal  loading  of  the  blade  and  the 
attachment. 

3.  The  gas  loads  on  the  blade  were  not  simulated  in  the 
test  rig. 

4.  The  difference  between  the  disk  and  blade  material 
moduli  (24  x  103  psi  for  the  disk  and  39  x  10®  psi 
for  the  blade)  was  not  simulated  in  the  photoelastic 
test  rig. 

A.*.i  of  the  abovft  would  have  contributed  to  the  differences 
observed  between  the  photoelastic  results  and  the  calculated 
s trasses.  The  ability  to  study  the  relative  effects  of  tol¬ 
erances  and  load  sharing  in  the  photoelastic  test  rig  made 
the  photcclasfcic-  testing  a  valuable  tool  in  establishing  the 
final  fir-tree  configuration. 
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3. 3. 1.3  Rotor  Disk 

During  the  preliminary  design  mechanical  studies,  calcula¬ 
tions  were  made  to  determine  the  temperature  distributions 
(Paragraph  2.8)  and  the  centrifugal  and  thermal  elastic 
stresses  (Paragraph  2.9)  of  the  rotor  disks,  it  was  apparent 
from  these  preliminary  analyses  that  the  magnitude  of  the 
thermal  stresses  presented  a  problem  that  required  solution. 
Studies  were  initiated  to  determine  a  cooling  scheme  which, 
without  increasing  the  cooling-air  flow  rate,  would  provide 
greater  disk  cooling  and  produce  a  more  favorable  tempera¬ 
ture  gradient  in  the  disk  to  reduce  the  thermal  stresses  to 
an  acceptable  level.  From  these  studies,  an  optimized 
cooling  scheme  was  determined  in  which  equal  portions  of 
the  cooling  air  were  impinged  upon  each  side  of  the  disk 
wall,  near  the  bore  of  the  disk,  to  flow  out  radially  along 
the  disk  sides  and  enter  the  hot-gas  stream  at  the  OD  of 
the  disk.  The  preliminary  estimates  for  the  airflow  rates 
were  determined  to  be  adequate — namely,  50  percent  for  the 
first  stage,  30  percent  for  the  second  stage,  and  20  percent 
for  the  third  stage,  with  the  total  representing  2.0  per¬ 
cent  of  the  engine  through-flow. 

With  completion  of  the  rotor  blade  and  the  rotor  disk/blade 
attachment  final  designs  and  with  the  cooling  scheme  opti¬ 
mized,  analyses  were  conducted  to  calculate  the  temperature 
distributions,  centrifugal  stresses,  and  the  combined  centrif¬ 
ugal  and  thermal  elastic  stress  of  the  first-stage  rotor 
disk. 

The  heat-transfer-coefficient  and  the  relative  gas- temperature 
distributions  around  the  blade  (Figure  28)  were  included  in 
the  final  thermal  analyses.  The  cooling-air-heat- transfer 
coefficients  were  calculated  with  use  of  the  estimated 
cooling-gas  temperatures.  A  coolant  inlet  temperature  of 
900°F  was  used  (this  value  was  based  on  a  5 . 0-pound- per- 
second  compressor  with  a  pressure  ratio  of  10:1,  and  an 
estimated  100  °F  temperature  rise  in  the  coolant  temperature 
from  the  time  it  leaves  the  compressor  discharge  until  it 
is  impinged  upon  the  sides  of  the  rotor  disk) .  An  existing 
computer  program  was  used  to  calculate  the  metal  temperatures, 
along  with  the  cooling-air  temperatures,  as  the  coolant 
flowed  radially  out  along  the  sides  of  the  disk. 

Thermal  studies  were  completed  for  the  rotor  with  and  without 
the  nickel  shim  stock  placed  between  the  mating  surfaces 
of  the  blade  and  disk  in  the  fir-tree  attachment.  As 
in  the  preliminary  studies  (Paragraph  2.8),  use  of  the 
shim  stock  was  determined  to  be  advantageous?  the  metal 
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temperatures  of  the  disk  in  the  attachment  area  were  reduced 
approximately  35°.  This  temperature  reduction  in  the  disk 
metal  at  the  attachment  area  will  allow  an  increase  of  approx¬ 
imately  5000  psi  in  stress  level  for  the  same  disk  life  with 
the  IN-100  disk  material.  In  addition,  use  of  the  shim  stock 
caused  a  slight  reduction  (1000  to  3000  psi)  in  the  peak 
thermal  stress  levels  in  the  disk.  Therefore,  the  shim  stock 
was  recommended  for  inclusion  in  the  final  design. 

In  the  temperature  analyses,  the  contact  resistance  across 
the  fir-tree  attachment  was  simulated  as  an  additional  length 
to  the  conduction  path  connecting  the  blade  to  the  disk. 

This  additional  conduction  length  was  estimated  (from  scant 
data  available  at  the  time)  as  0.09  inch.  In  terms  of  thermal 
resistance  across  the  fir-tree  attachment  contact  surfaces, 
this  value  indicated  that  the  resistance  was  quite  small  and 
would  have  very  little  effect  on  the  thermal  gradients  be¬ 
tween  the  blade  and  the  disk.  There  was  concern  that  this 
resistance  might  be  of  much  larger  magnitude;  therefore,  a 
test  program  was  prepared  to  determine  the  thermal  conduction 
across  the  interface  of  the  fir-tree  attachment  between  the 
blade  and  the  disk. 

The  thermal  conduction  test  program,  reported  in  Appendix  IV 
of  this  document,  consisted  of  two  major  parts s 

1.  Design  and  fabrication  of  the  thermal  conduction 
test  rig  and  test  specimens 

2.  Implementation  of  the  test  and  analysis  of  the 
data 

Results  from  the  thermal-conduction  testing  indicated  that 
the  contact  resistance,  instead  of  being  larger,  was  approx¬ 
imately  50  percent  smaller  than  that  used  in  rotor  thermal 
analyses.  Since  the  contact  resistance  was  so  small,  reducing 
it  further  would  have  little  or  no  effect  on  the  calculated 
disk  temperatures  and  cooling  requirements,  so  no  additional 
disk  thermal  analyses  were  conducted. 

The  calculated  temperature  distributions  of  the  rotor  with 
and  without  use  of  the  shim  stock  are  shown  in  Figure  33. 

The  calculated  disk  and  blade  centrifugal  stresses  and  the 
combined  centrifugal  and  thermal  elastic  stresses  of  the 
rotor,  without  use  of  the  shim  stock,  are  shown  on  Figures  34 
and  35,  respectively.  These  stress  data  were  selected  to 
illustrate  that  the  disk  stresses  are  within  an  acceptable 
range  even  without  the  benefits  gained  by  use  of  the  shim 
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Figure  33.  Temperature  Distributions  of  Rotor 

With  and  Without  Shim  Stock  in  Place 
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Figure  34.  Calculated  Centrifugal  Stresses  of 
Rotor  Disk  Without  Shim  in  Place, 


78 


1.  IMPINGEMENT  COOLANT  FLOW  SPLIT  EQUALLY  TO  EACH  SIDE 
OF  THE  DISK  WITH  TOTAL  COOLANT  FLOW  EQUAL  TO  1- 
PERCENT  OF  ENGINE  THROUGH -FLOW 

2.  TURIINE  SPEED  =  3R.000  RPM 

3.  FIRST-STAGE  TURIINE  INLET  TEMPERATURE  =  2400°F 

4.  COOLANT  INLET  TEMPERATURE  =  900°F 

5.  POINTS  MARKED0  REFER  TO  DESIGN  FILLET  STRESSES  IN  DISK 
PORTION  OF  HADE- DISK  ATTACHMENT 

6.  POINTS  MARKED  O  REFER  TO  DESIGN  FILLET  STRESSES  IN 
•LADE  PORTION  OF  ILAOE-DISK  A  ATTACHMENT 

7.  COOLANT  FLOW  IMPINGES  AT  RADIUS  =  1.400  INCH 
I.  ENGINE  TOTAL  THROUGH-FLOW  =  S  LIS  PER  SEC. 

9.  DISK  A  AVERAGE  TANGENTIAL  STRESS  =  54.  500  PSI 
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Figure  35.  Calculated  Combined  Centrifugal  and 
Thermal  Elastic  Stresses  of  Rotor 
Disk  Without  Shim  in  Place. 
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stock.  The  combined  centrifugal  and  thermal  stresses  in  the 
disk  bore  (Figure  35)  may  result  in  some  yielding  in  that 
area.  However,  successful  experience  with  similarly  shaped 
disks  operating  at  the  same  stress  levels  and  in  a  similar 
temperature  environment  provides  assurance  that  the  magnitude 
of  the  predicted  bore  stresses  will  not  result  in  low-cycle- 
fatigue  problems. 

3.3.2  Stator — Mechanical  Design  Studies 

Detailed  stress  analyses  were  conducted  on  *\e  final  aero¬ 
dynamic  vane  shape.  These  studies  included  calculations  of 
the  maximum  transient  thermal  stresses  that  would  occur  in 
the  vane  during  an  engine  start  and  also  the  peak  bending 
stresses  that  would  occur  in  the  vane  from  the  steady -state 
air  loads.  The  results  of  the&e  studies  were  to  be  used  to 
establish  the  strength  goal's  of  the  beryllide  material  (Com¬ 
position  1)  to  be  developed  for  stator  -application. 

v 

The  same  analytical  procedures,  assumptions,  and  computer 
programs  that  were  used  for  the  rotor  blade  thermal  studies 
(Paragraph  3.3. 1.1)  were  used  for  the  temperature  and  thermal- 
stress  analyses  of  the  stator  vane.  The  beryllide  Composi¬ 
tion  1  physical  properties  that  were  estimated  for  these 
studies  were  shown  in  Table  VIII.  While  thermal-stress 
analyses  were  performed  on  the  stator  vane  designs  under  con¬ 
sideration  for  four  different  engine  transient  temperature 
profiles,  the  data  presented  in  this  report  were  based  on  the 
same  assumed  TIT-variation-with-time  temperature  profile  that 
was  used  for  the  rotor  blade  analyses.  This  temperature 
profile  (Figure  27)  was  considered  to  be  the  most  realistic 
estimate  of  the  transient  temperature  conditions  that  would 
be  encountered  during  a  normal  start  of  an  engine  of  the 
design,  configuration,  and  size  assumed  for  this  program. 
Calculations  of  the  peak  bending  stresses  in  the  vane  due  to 
the  air  loads  were  based  on  the  static  pressure  distribution 
around  the  vane  calculated  from  the  vane  velocity  distribu¬ 
tions  established  in  the  aerodynamic  final  design  studies 
(Figure  23) . 

Since  the  stator  vane  has  a  fairly  constant  cross  section  with 
radius,  the  data  presented  for  the  stator  studies  are  for  the 
cross  section  located  at  8  =  3.5694  and  are  considered  to  be 
representative  values  for  all  cross  sections.  Subsequent 
paragraphs  describe  the  pertinent  mechanical  studies  that 
were  conducted  for  the,  stator. 
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3 . 3 . 2 . 1  Two-Dimensional  Thermal  and  Thermal-Stress  Transient 
Analyses  of  the  Stator  Vane 

These  analyses  comprised  several  phases,  which  are  described 
In  chronological  order  on  the  following  pages.  The  initial 
study  was  undertaken  to  calculate  the  heat-transfer  coeffi¬ 
cient  and  relative  gas- temperature  distributions  around  the 
vane  at  steady-state  conditions.  The  calculations  were  based 
on  the  gas  velocity  distributions  (Figure  23)  and  cycle  inlet 
conditions.  These  cycle  conditions,  along  with  the  resulting 
heat-transfer  coefficient  and  relative  gas- temperature  dis¬ 
tributions,  are  shown  in  Figure  36. 

The  first  mechanical  configuration  of  the  aerodynamic  vane 
design  to  be  evaluated  was  a  solid  vane.  Results  of  this  two- 
dimensional  thermal  and  thermal-stress  study  are  shown  in 
Fioure  37.  The  peak  thermal  stress  occurred  as  a  compressive 
stress  near  the  leading  edge  of  the  vane.  The  thermal  stress 
and  temperature  distributions  in  the  solid  vane  at  the  ap¬ 
proximate  time  of  maximum  stress  are  shown  in  Figure  38. 

The  maximum  thermal  stresses  for  the  solid  vane  are  of  the 
same  order  of  magnitude  as  those  calculated  for  the  hollow 
first-stage  blade  discussed  in  Section  3. 3. 1.1  and,  thus,  were 
considered  acceptable.  During  the  analysis  it  was  observed, 
however,  that  the  magnitude  of  the  maximum  thermal  stress 
levels  in  the  stator  vane  were  due,  in  part,  to  the  large 
variation  in  gas  velocities  around  the  vane  and  the  corres¬ 
ponding  heat-trans f er  coefficient  distribution,  it  was  ap¬ 
parent  that  by  discretely  locating  a  radial  contoured  hole  in 
the  vane,  these  maximum  stresses  could  be  reduced.  Since  in 
the  absence  of  beryllide  material  properties,  the  goal  was  to 
minimize  the  stress  levels  occurring  in  the  beryllide  parts, 
analysis  was  undertaken  to  determine  the  location  and  contour 
of  the  radial  hole.  Calculations  were  performed  to  determine 
a  ratio  of  pressure- side-wall  thickness  to  suction-side-wall 
thickness  that  would  serve  to  maintain  a  minimum  temperature 
difference  between  the  two  vane  walls  during  transient  condi¬ 
tions.  The  configuration  established  was  an  offset-hollow 
vane  with  wall  thicknesses  of  0.02Q  inch  for  the  cool  wall 
(pressure  side)  and  0.50  inch  for  the  hot  wall  (suction 
side) . 


A  two-dimensional  thermal  and  thermal- stress  study  was  then 
performed  on  the  offset-hollow  vane.  Results  of  this  study 
are  shown  in  Figure  39.  The  peak  compressive  stress  occurred 
on  the  pressure  side  of  the  vane  just  aft  of  the  radial  hole 
centerline,  and  the  peak  tensile  stress  occurred  near  the 
leading  edge  of  the  vane.  The  thermal  stress  and  temperature 
distributions  in  the  offset-hollow  vane  at  the  approximate 
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Figure  36.  Steady-State  Heat-Transfer  Coefficient  and 
Relative  Gas -Temperature  Distributions 
Around  the  Stator  Vane  at  Section, 

R  =  3.5694  inches. 
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Figure  37.  Transient  Thermal  Study  of  Solid 

Stator  Vane  at  Section,  R  =  3.5694  Inches. 
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Figure  39.  Transient  Thermal  Study  of 

Offset-Hollow  Stator  Vane  at 
Section,  R  =  3.5694  Inches. 
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Figure  40.  Stress  and  Temperature  Distributions 
in  Offset-Hollow  Stator  Vane  at 
Section ,  R  *  3,5694  Inches. 


time  of  maximum  stress  are  shown  in  Figure  40.  The  offset- 
hollow  vane  configuration  resulted  in  an  approximate  30  per¬ 
cent  reduction  of  the  peak  thermal  stresses  that  were  cal¬ 
culated  for  the  solid  vane. 

3. 3. 2. 2  Stator-Vane  Bending  Stress  Studies 

Studies  were  made  to  determine  the  peak  steady-state  bending 
stresses  in  both  the  solid  and  offset-hollow  vanes  that  result 
from  the  static  pressure  distribution  around  the  vane.  The 
calculated  peak  bending  stresses  were  11,000  psi  for  the  solid 
vane  and  14,000  psi  for  the  offset-hollow  vane.  For  both  vane 
configurations,  the  peak  bending  stress  was  a  tensile  stress 
located  <at  the  leadihg  edge  of  the  vane  where  it  intersects 
the  outer^cifcumf er erjce  of  the  nozzle  passage.  Shear  stresses 
in  both  vane  configurations  due  to  static  pressure  loading 
were  calculated  to  be  less  them  1000  psi.  Figure  41  shows 
the  location -And  position  of  the  principal  axis  of  the  vane 
for  both  vane  configurations  used  in  this  study. 

These  steady-state  bending  stresses  in  the  vane  were  used  to 
establish  the  stress-rupture  requirement  for  the  beryllide 
vane  material. 

3. 3. 2. 3  Stator- Vane  Attachment  Studies 


At  this  time,  the  anticipated  mechanical  configuration  for 
the  stator  assembly  was  a  series  of  one-piece  stator  vanes 
attached  to  an  outer  ring  to  form  the  assembly.  Two  mechan¬ 
ical  stator-vane  attachment  configurations  were  designed  for 
evaluation.  These  configurations  were  identified  as  A  and 
B,  and  are  shown  in  Figure  42. 

The  static  pressure  drop  across  the  stator  vane  was  calculated 
to  be  60  psi.  This  degree  of  pressure  differential  imposed 
considerable  vane  loading  on  the  attachment  and  also  presented 
a  serious  sealing  problem  in  the  vane  attachment.  Therefore, 
both  attachment  configurations  were  designed  to  ensure  ad¬ 
equate  sealing  and  to  keep  the  maximum  stresses  in  the  attach¬ 
ment,  resulting  from  the  static  pressure  loading  on  the  vane, 
to  less  than  the  steady-state  bending  stresses  occurring  in 
the  vane. 

Both  attachment  configurations  were  designed  for  use  of  the 
labyrinth  provided  by  lap  joints  for  scaling  between  the 
individual  vanes.  Also,  piston  ring  arrangements  similar  to 
those  used  on  the  Model  TPE331  Turboprop  Engines  were  designed 
for  the  stator  assembly  ID  to  allow  for  small  axial  movements 
and  still  provide  a  good  method  of  sealing.  The  piston  rings 
would  provide  an  elastic  member  to  control  the  deformations 
of  the  individual  vanes. 
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Figure  41.  Locations  and  Positions  of 

Principal  Axis  on  Offset-Hollow 
and  Solid  Stator  Vanes  at 
Section,  R  =  3.5694  Inches. 


Two  separate  methods  were  designed  to  provide  sealing  between 
the  stator  shroud  and  the  attachment  ring.  Configuration  A 
was  designed  for  use  of  a  flat  sealing  ring  clamped  in  place* 
Configuration  B  again  employed  the  labyrinth  of  a  lap  joint. 

The  plan  at  this  time  was  to  evaluate  both  configurations 
experimentally  during  the  cascade  test  program  scheduled  for 
the  stator  component.  Based  on  the  results  of  the  cascade 
testing,  one  configuration  would  then  be  selected  for  the 
final  design. 

3. 3.2.4  Three-Dimensional  Thermal  Analyses  of  the  Offset- 
Hollow  Stator  Vane 

With  mechanical  configurations  of  the  stator-vane  attachment 
established,  three-dimensional  thermal  analyses  on  the  offset- 
hollow  vane  with  attachment  Configuration  B  were  initiated. 
Both  steady- state  and  transient  analyses  were  included. 

The  same  computer  programs  and  analytical  procedures  used 
for  the  two-dimensional  analyses  were  used  for  the  three- 
dimensional  studies  plus  the  following  inclusions  and  assump¬ 
tions: 

1.  Radial  heat  conduction  in  the  vane  was  included 
in  the  study. 

2.  The  effects  of  compressor  discharge  air  scrubbing 
in  the  nozzle  vanes  supporting  ring  were  included 
in  the  study. 

3.  For  the  steady-state  thermal  study,  all  the  heat 
flowing  into  the  vane  was  assumed  to  flow  out 
through  the  vane  attachment  to  the  supporting  ring 
located  on  the  OD  of  the  stator  assembly. 

4.  For  the  transient  thermal  study,  the  assumed  initial 
temperature  of  all  metal  components  was  80°F.  The 
assumed  compressor  discharge  temperature  variation 
with  time  is  shown  in  Figure  43.  The  TIT  variation 
with  time  used  in  the  preceding  thermal  studies  was 
again  used  for  this  study  (Figure  27) . 

The  calculated  steady-state  temperature  distribution  for  the 
offset -hollow  vane  and  attachment  Configuration  B  is  shown  in 
Figure  44.  The  steady-state  thermal  study  revealed  that  most 
of  the  nozzle  vane  was  at  or  near  the  relative  gas  tempera¬ 
ture.  It  was  found  that  a  rather  large  thermal  gradient 
existed  in  the  stem  of  the  attachment  portion  of  the  vane.  It 
was  determined  that  this  thermal  gradient  could  present  an 


90 


Figure  43.  Assumed 


Figure  44.  Three-Dimensional,  Steady-State 
Temperature  Distribution  in 
Offset-Hollow  Vane  and  Attachment 
Configuration  B. 


unacceptable  stress  problem  with  the  one-piece  vane  design 
due  to  the  bending  imposed  on  the  platform  by  the  stem  or 
dovetail  and  the  restraint  provided  by  the  vane  itself.  At 
this  time,  studies  were  initiated  to  design  a  multiple-piece 
vane  (described  in  3. 3. 2. 5)  that  would  help  to  relieve  these 
stresses  as  well  as  ease  the  fabrication  problems  encountered 
in  fabricating  a  one-piece  vane. 

Results  of  the  transient  thermal  study  were  used  to  evaluate 
the  accuracy  of  the  two-dimensional  thermal  study  reported 
in  3. 3. 2.1.  Figure  45  shows  a  comparison  of  the  calculated 
temperatures  from  the  two-dimensional  study  assuming  no 
radial  heat  conduction  and  from  this  three-dimensional  study 
assuming  radial  heat  conduction.  The  temperatures  shown  in 
Figure  45  are  for  vane  locations  where  the  maximum  calculated 
thermal  stresses  occurred  in  the  offset-hollow  vane.  There 
is  good  agreement  between  the  two  analyses  up  to  4.0  seconds, 
which  is  well  beyond  2.5  seconds — the  approximate  time  of 
peak  thermal  stress.  Therefore,  it  was  substantiated  that 
the  two-dimensional  thermal  studie,$.'*accurately  predicted  the 
maximum  transient  temperature  gradients. 

3. 3. 2. 5  Multiple-Piece  Stator-Vane  Studies 

Because  of  difficulties  encountered  in  fabricating  a  single 
or  one-piece  stator  vane  (platforms  and  attachment  integral 
with  the  airfoil  portion) ,  it  was  decided  to  fabricate  each 
vane  insert  as  three  pieces.  The  insert  was  to  include  the 
platform  and  attachment  portion  on  the  outer  diameter  as 
one  piece,  the  airfoil  as  one  piece,  and  the  platform  and 
attachment  portion  located  on  the  inner  diameter  as  the  third 
piece.  Following  a  review  of  the  loading  on  the  vane,  it 
was  determined  that  this  change  would  also  improve  the  stress 
situation  in  the  vane.  This  improvement  was  brought  about 
by  the  method  now  required  to  support  the  three-piece  vane. 

Since  it  would  now  be  an  assembly  of  three  pieces,  provision 
for  support  of  the  vane  both  on  the  inner  diameter 
and  the  outer  diameter  was  required.  The  previous  one-piece 
stator  vane  was  essentially  cantilevered  from  a  support  on 
its  outer  diameter.  In  the  three-piece  configuration,  the 
vane  was  simply  supported  at  each  end,  which  reduced  the 
steady-state  bending  stresses  to  approximately  one- fourth 
of  their  previous  value  in  the  one-piece  design.  The  steady- 
state  bending  stresses  in  the  airfoil  for  the  three-piece 
vane  were  calculated  to  be  2500  psi  in  tho  solid  vane  and 
3250  psi  in  the  hollow  vane. 
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NOTES: 

1.  CURVES  REPRESENT  DATA  FROM  TWO- 
DIMENSIONAL  ANALYSIS 

2.  DATA  POINTS  REPRESENT  RESULTS 
FROM  THREE-DIMENSIONAL  ANALYSIS 

3.  POINTS  l  AND  2  ARE  POINTS  OF  MAXIMUM  STRESSES 


TIME,  SEC. 


Figure  45.  Comparison  of  Two-Dimensional  and 
Three-Dimensional  Calculated 
Temperatures  for  Offset-Hollow  Vane 
at  Points  of  Maximum  Stresses. 
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Based  on  the  results  of  the  three-dimensional  heat  study 
(Section  3. 3. 2. 4),  no  change  was  expected  in  the  transient 
thermal  stress  results  presented  in  Section  3. 3. 2.1. 

3. 3. 2. 6  Summation  of  the  Stator  Mechanical  Studies  and 
Cascade  Test  Evaluation 


When  the  material- property  data  became  available  for  beryllide 
Composition  1,  it  wa3  found  to  be  inadequate  to  accommodate 
the  stresses  occurring  in  the  three-piece  vane.  As  reported 
in  Volume  I,  this  was  verified  by  the  cascade  testing  where 
the  beryllide  test  vanes  failed  in  a  relatively  short  period 
of  time.  These  failures  were  primarily  caused  by  insuffi¬ 
cient  low-temperature  ductility  in  the  beryllide  to  tolerate 
the  peak  thermal  stresses  occurring  during  starts  and  calcu¬ 
lated  in  the  mechanical  studies.  In  addition,  the  material- 
property  test  results  showed  that  the  beryllide  had  inadequate 
stress-rupture  strength  to  satisfy  the  steady-state  stress 
requirements  of  the  vane. 

Although  the  vane  steady-state  stresses  had  been  compared  to 
the  stress-rupture  strength  of  the  beryllide  for  200  hours 
at  2400°F,  a  more  realistic  approach  is  to  compare  the 
stresses  to  the  100-hour  0.2-percent  creep  strength.  This 
comparison  was  made  by  estimating  the  0.2-percent  creep 
strength  to  be  equal  to  80  percent  of  the  stress-rupture 
strength.  Thus,  the  stress-rupture  strength  required  for 
the  three-piece  solid-vane  configuration  becomes  3130  psi 
(the  calculated  stress  divided  by  0.8)  for  100  hours.  The 
temperature  used  in  this  evaluation  includes  a  ±200°F  allow¬ 
ance  for  the  circumferential  spread  in  combustor  discharge 
temperature.  This  gives  a  maximum  expected  vane  metal  tem¬ 
perature  of  2600°F.  With  this  temperature  and  the  stress- 
rupture  strength  presented  in  Volume  I,  the  0.2-percent 
creep  life  of  the  Composition  I  solid  vane  is  estimated  to 
be  1.25  hours. 
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4.  CONCLUSIONS 


When  the  material-property  tests  (reported  in  Volume  I)  were 
completed#  the  resulting  data  were  compared  with  the  turbine 
design  requirements  calculated  in  the  final  design  mechanical 
studies.  The  following  conclusions  were  made: 

1.  The  beryllide  material  had  insufficient  ductility 
for  use  in  either  the  stator  or  the  rotor  components. 
Measurable  ductility  in  the  beryllides  was  not 
achieved  below  2000°F. 

2.  The  beryllide  material  had  insufficient  stress- 
rupture  strength  for  use  in  the  rotor  component. 

While  Composition  2  displayed  adequate  stress- 
rupture  strength  for  the  stator  design  (200  hours  at 
2500  psi  and  2400°F),  it  did  not  approach  the  rotor 
design  requirement. 

Ir  view  of  these  findings,  it  was  decided  that  without  addi¬ 
tional  research  to  further  develop  a  high-temperature  turbine 
material#  it  vas  not  expedient  to  continue  design  studies  for 
the  uncooled  turbine. 
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RECOMMENDATIONS 


The  objective  of  this  activity — to  demonstrate  the  feasibil¬ 
ity  of  an  uncooled- turbine  design — was  not  realized.  Success 
of  the  program  was  precluded  by  the  lack  of  a  suitable  tur¬ 
bine  material.  Even  so,  the  need  for  such  a  turbine  becomes 
increasingly  greater,  and  the  successful  design  of  a  small, 
high- temperature  uncooled  turbine  should  remain  an  ultimate 
goal  advanced  turbine  technology. 

While  it  is  advocated  that  no  further  turbine  design  effort 
be  expended  at  this  time,  it  is  recommended  that  research  o 
high- temperature  materials  be  continued  in  an  effort  to  de¬ 
velop  a  material  that  would  facilitate  the  success  of  an 
'  uncooled-turbine  design. 
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APPENDIX  I 

GENERATION  OF  PRELIMINARY  BLADE  PROFILES 


Figure  46  shows  a  representative  blade  profile,  along  with 
the  definition  of  the  geometry  used.  The  axial  chord  CX  and 
blade  spacings  are  obtained  by  Zweifel's  solidity  require¬ 
ment  (reference  on  Page  19) : 


CX 

s 


2  cos  83  sin  (P2  -  63> 


(24) 


where,  for  this  study,  Y  =  0.8  (Zweifel's  recommended  opti¬ 
mum  value) . 

From  this  ratio. 


where  r^  =  radius  of  blade  section 

Z  =  number  of  blades 


(25) 


The  blade  profile  layout  (see  Figure  46)  is  obtained  by  lo¬ 
cating  the  center  of  a  single  arc  suction-surface  radius  r 

\  ® 

so  that  the  longest  possible  arc  joins  a  tangent  vith  blade 
angle  &3  at  the  exit  and  a  tangent  with  angle  B2  at  the  blade 

leading  edge.  This  suction-surface  radius  is  given  by 


CX  +  s  sin  63  cos  f?3 


(26) 


The  blade  pressure  surface  is  also  a  single  circular  arc  ob¬ 
tained  by  locating  a  circle  through  two  points  and  tangent 
to  a  trailing-edge  surface  t  thick.  The  first  point  estab¬ 
lishing  the  pressure  surface  is  obtained  by  laying  out  a 
leading-edge  thickness  dimension  t  from  the  leading-edge 

intersection  of  the  suction-surface  arc.  The  second  point  on 
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the  pressure  surface  is  obtained  by  laying  out  a  camber 
thickness  TC  on  a  suction-surface  radius  at  a  point  halfway 
through  the  suet ion- surface  turning  angle  (82  -  83). 

A  computer  program  calculates  these  blade  layout  dimensions 
and  the  blade  areas  and  blade  moments,  and  integrates  them  to 
give  blade  volumes  and  radial  centroids  from  input  data  82, 

8~,  Z,  *■  ,  TC/CX,  t  ,  and  radius  of  the  blade  section  r.. 

For  the  turbines  studied  in  this  program,  t  =  0.070  and 

a 

t  =  0.030.  These  values  were  judged  to  be  reasonable  from 
a  fabrication  point  of  view. 

Since  the  blade  areas  are  a  function  of  certain  input  param¬ 
eters  such  as  number  of  blades  or  assumed  thickness  parame¬ 
ters  t  ,  TC/CX,  and  t  ,  a  study  of  blade-area  variation  with 

3  G 

two  of  these  parameters  was  made  for  the  third  stage  of  Tur¬ 
bine  Design  A-338-5. 

Figure  47  shows  blade-area  variation  with  radius  for  35,  40, 
and  45  blades.  Figure  48  shows  the  effect  of  camber  thick¬ 
ness  -on  blade  area.  These  two  figures  were  prepared  for  use 
in  the  stress  analysis  to  show  sensitivity  to  these  assump¬ 
tions  in  the  blade  profile  layout. 
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Figure  47.  Turbine  Blade  Area  Distribution  Turbine  Design 
A-338-5,  Third  Stage,  Showing  Comparison 
Between  Different  Numbers  of  Blades. 
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APPENDIX  II 

METHOD  FOR  ESTIMATING  TURBINE 
AERODYNAMIC  EFFICIENCY 


In  general,  turbine-blade  losses  can  be  categorized  as  two- 
dimensional  profile  losses  and  three-dimensional  losses.  The 
two-dimensional  losses  consist  of  surface  friction  losses  due 
to  boundary  layers  and  downstream  mixing  losses  due  to 
boundary- layer  and  trailing-edge  wakes.  The  three-dimensional 
losses  consist  of  end-wall  boundary  layers,  secondary  flow 
losses,  and  tip  clearance  losses.  Two-dimensional  profile 
losses  can  be  predicted  with  reasonable  accuracy  as  long  as 
the  wake  region  is  so  thin  that  it  has  no  effect  on  the  pro¬ 
file  velocity  distribution.  Three-dimensional  profile  losses 
present  a  more  difficult  problem,  and  experimental  data  are 
required. 

For  a  wide  variety  of  turbine  blading,  Stewart,  Whitney,  and 
Wong*  have  shown  that  the  trailing-edge  momentum  thickness 
can  be  correlated  with  total  blade  surface  diffusion  by  the 
following  equations 


0.003 
1-1. 4D. 


x  10‘ 
Re, 


0.2 


(27) 


where  8  =  total  trailing-edge  momentum  thickness 

l  =  mean  camber-line  length  of  blade  section 

Re^  =  blade  row  Reynolds  number  based  oh  camber¬ 
line  length 

D^  =  total  surface  diffusion  parameter 


*  Stewart,  Warner  L. ,  Whitney,  Warren  J. ,  and  Wong,  Robert  Y. , 
A  STUDY  OF  BOUNDARY  LAYER  CHARACTERISTICS  OF  TURBOMACHINE 
BLADE  ROWS  AND  THEIR  RELATION  TO  OVERALL  BLADE  LOSS,  ASME 
Paper  No.  59-A-23,  1959. 
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The  total  diffusion  parameter  is  the  sum  of  the  suction  sur¬ 
face  and  pressure  surface  diffusion  parameters  defined  by: 


D 


D 


s 


=  1 


+  D 

P 

V 

exit 

V 

max 

Vmin 

v 

vmlet 


where 


=  total  diffusion  parameter 
Dg  =  suction-side  diffusion  parameter 

Dp  =  pressure-side  diffusion  parameter 

Vexit  =  averaSe  velocity 

V  =  maximum  velocity  on  the  suction  side 
max 

V  .  =  minimum  velocity  on  the  pressure  side 


These  parameters  are  related  to  the  blade  loading  diagram  as 
shown  in  Figure  49. 

In  addition,  Zweifel  (reference  on  Page  19)  has  shown  that 
the  optimum  blade  section  solidity  can  be  related  to  the 
Dlaae  geometry  by  tne  following  expression: 


CX 

s 


2  cos  3^  sinA3 
exit 

_  3. 

cos  inlet 


(28) 


where 


?  =  aerodynamic  loading  coefficient 


SURFACE  VELOCITY  V 
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Figure  49.  Blade  Aerodynamic  Loading  Diagrams 
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The  aerodynamic  loading  coefficient  Y  is  a  complex  function 
that  is  related  to  the  profile  suction  and  pressure-side 
velocities.  If  the  idealized  suction  surface  velocity  were 
constant  and  equal  to  the  exit  velocity  and  the  pressure  sur¬ 
face  velocity  were  zero,  the  loading  coefficient  would  be 
unity.  The  parameter  Y  varies  from  0.65  to  0.80  for  highly 
loaded  blades. 

With  the  assumption  of  an  idealized  rectangular  loading  dia¬ 
gram  having  zero  suction-surface  diffusion  (D  =0),  the 

s 

parameter  Y  can  be  related  to  the  total  diffusion  parameter 
Dt  by  the  following  expression  (see  Figure  49) : 


By  combining  the  above  equations,  the  two-dimensional  surface 
friction  losses  can  be  related  to  the  blade  geometry  and  the 
requirements  of  the  vector  diagram.  With  use  of  the  above 
information,  an  estimate  of  the  wake  mixing  losses  can  be 
made  according  to  the  Stewart  method.*  This  method  can  be 
refined  if  detailed  blade  sections  are  available,  since  the 
potential  flow  velocity  distribution  about  the  blade  profile 
can  be  calculated  and  the  results  used  to  determine  the 
boundary- layer  growth  on  the  blade  surfaces. 

Integration  of  the  estimated  two-dimensional  blade  section 
losses  over  the  entire  blade  height  and  correction  of  the 
results  for  end-wall  boundary  layers  provide  an  estimated 
loss  for  the  blade  row.  The  end-wall  boundary- layer  cor¬ 
rection  is  described  by  Stewart,  Whitney,  and  Wong 
(reference  on  Page  103) .  Combining  stator  and  rotor  loss 
estimates  results  in  a  predicted  aerodynamic  stage 
efficiency,  a- 


♦Stewart,  Warner  L.,  ANALYSIS  OF  TWO-DIMENSIONAL  COMPRESSIBLE 
FLOW  LOSS  CHARACTERISTICS  DOWNSTREAM  OF  TURBOMACHINE  BLADE 
ROWS  IN  TERMS  OF  BASIC  BOUNDARY  LAYER  CHARACTERISTICS,  NACA 
TN  3515,  July  1955. 
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Tip-clearance  losses  in  the  rotor  are  correlated  by  the  fol¬ 
lowing  equation: 


"  ^3D,  A 


rtk5 

1  -  *  h 

m 


(30) 


where 


^3D  =  stage  efficiency  for  rotor- tip  clearance 

equal  to  c 

tj- _  =  predicted  three-dimensional  aerodynamic 

A  efficiency  at  zero  clearance 


V  rm 


=  rotcr-tip  and  mean  radius,  respectively 


c  =  rotor-tip  clearance 

h  =  rotor-blade  height 

k  =  constant,  depending  on  tip  shroud  geometry 

Figure  50  shows  a  typical  loss  breakdown  for  a  single-stage 

turbine  for  varying  rotor  clearance  and  tip  shroud  geometry. 
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TURBINE  EFFICIENCY 


Figure  50.  Attainable  Turbine  Efficiency. 
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GENERAL 


A  photoG Id s t ic  test  program  was  implement^'.,  t x  c-xper  1. -c: a  !  1  . 
evaluate  the  stress  distribution  of  the  rotor  Jins  blalc 
fir-tree  attachment  configurations  that  evolved  iroro  the 
mechanical  studies  reported  in  Paragraph  3. 3. 1.2. 2.  The 
attachment  configurations  were  to  be  tested  over  a  wide  range 
of  loads,  simulating  operating  conditions,  while  observing 
the  corresponding  stress  distributions.  Information  gained 
from  these  tests  was  to  be  used  in  studies  to  optimize  the 
rotor  final  mechanical  design  with  respect  to  stress  concen¬ 
tration  in  the  attachment  area. 

The  model  material  selected  for  the  photoelastic  test  evalu¬ 
ation  was  CR-39,  an  ailyl  diglycoi  carbonate  resin.  For  the 
results  of  the  photoelastic  testing  to  be  applicable  to  the 
actual  turbine  prototype,  it  was  necessary  that  certain 
conditions  of  structural  similarity  be  satisfied.  Conditions 
of  similitude  were  established  through  dimensional  analyses, 
and  the  model  loads  determined  were  used  during  the  test 
evaluations. 

SIMILITUDE  CONDITIONS 

The  method  for  establishing  similitude  conditions  between 
models  and  prototypes  by  dimensional  analyses  was  used  for 
the  photoelastic  test  model. 

With  the  assumption  that  both  the  model  and  the  prototype 
behave  elastically,  and  with  consideration  that  for  the 
general  case  the  load  distribution  may  be  affected  by  the 
elastic  deformation,  it  was  found  that  the  conditions  of 
similitude  were  governed  by  the  relationship 


(31) 


where  the  subscripts  m  and  p  refer  to  the  model  and  the  proto¬ 
type,  respectively;  thus. 
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y-  =  scaling  factor  in  the  plane  of  the  wheel 
P 

t 

—  -  scaling  factor  in  the  axial  direction 
~P 

For  convenience,  the  scaling  factors  were  made  5:1  in  the 
plane  of  the  wheel  and  1:4  axially.  The  prototype  (centrif¬ 
ugal  force)  load  at  operating  speed  was  calculated  to  be 
2500  pounds.  Available  data  indicated  a  Young's  modulus  of 
between  10  x  106  and  30  x  10e  psi  for  the  beryllide  proto¬ 
type  material  and  0.3  x  106  psi  for  the  CR-39  model  material. 
The  requisite  model  loads  were  then  calculated  to  be  between 
31  and  94  pounds,  depending  on  the  elastic  modulus  of  the 
beryllide. 

If  the  structure  were  sufficiently  stiff,  the  stress  distri¬ 
bution  would  be  unaffected  by  the  deformations,  and  the 
preceding  restriction  on  model  loads  would  be  unnecessary. 

CALIBRATION  TEST 

The  stress-optical  characteristics  of  the  CR-39  model  mater¬ 
ial  were  determined  by  a  calibration  test  of  a  representative 
tensile  specimen.  The  fringe  value  for  the  material  thick¬ 
ness  of  l/4  inch  was  found  to  be  340  psi  per  fringe  for  the 
monochromatic  light  transmitted  through  the  filter  used. 

This  was  in  agreement  with  the  manufacturer's  published 
value  of  336  psi  per  fringe  for  light,  also  presumably  in 
the  blue-green  portion  of  the  spectrum. 

The  value  obtained  in  the  calibration  test  was  used  in  the 
conversion  of  fringe  value  to  stress  value  in  the  actual 
model  test. 

TEST  DESCRIPTION 

The  model  was  made  five  times  actual  size  and  was  patterned 
after  Figure  51.  Three  blade  sections  were  used  to  simulate 
conditions  in  the  vicinity  of  the  center  blade  on  which  all 


measurement s  were  sad-  . 


t  *  x  > 


t  * 


The  loading  trame  with  the  model  ir.  place  ir  ?h»-  polariscope 
field  is  illustrated  in  ligure  53. 

Load  was  applied  by  adding  water  to  the  initially  counter¬ 
balanced  tank  shown.  This  was  transmitted  to  the  model 
through  the  lever-fulcrum  arrangement  and  a  set  of  calibrated 
force  rings  that  were  instrumented  to  give  readings  of  the 
loads  on  each  blade.  I’he  force  rings  had  a  capacity  of 
200  pounds  each  and,  with  the  attendant  instrumentation, 
were  capable  of  giving  load  readings  to  within  ±0.5  pound 
over  the  entire  load  range.  A  fourth  force  ring  of  700-pound 
capacity  served  as  an  added  check  on  the  total  load.  The 
force  rings  were  threaded,  which  permitted  equalization  of 
load  on  each  of  the  blades. 

Because  of  the  low  fringe  orders  developed  in  the  model  at 
the  load  range  of  interest,  attempts  were  made  to  use  a 
quartz  wedge  compensator  to  determine  fractional  fringe 
orders.  However,  the  high-stress  gradient  in  the  region 
of  maximum  stress  resulted  in  a  loss  of  distinctiveness  in 
the  fringes.  As  a  result,  this  approach  was  abandoned. 

In  lieu  of  this,  the  polariscope  was  set  for  light-field, 
circularly  polarized  light,  and  the  blade  load  was  recorded 
when  the  l/2-  and  3/2-order  fringe  appeared  at  the  point 
of  maximum  stress.  Similarly,  the  loads ‘were  recorded  on 
appearance  of  the  first-  and  second-order  fringes  with  the 
use  of  a  dark-field  polariscope.  Such  load  readings  were 
taken  over  a  number  of  load  and  unload  cycles  of  the  specimen. 
The  time  between  cycles  (i.e. ,  recovery  period)  was  generally 
greater  than  1  hour. 

Finally,  as  an  additional  experiment,  the  model  was  loaded 
to  arbitrary  values  in  the  desired  load  range,  and  the 
fringe  values  were  determined  at  points  of  maximum  stress 
to  an  accuracy  of  approximately  1/25 th  of  a  fringe. 

PHOTOELASTIC  RIG  TESTING 


Upon  application  of  the  load  to  the  model,  the  lower  lobes 
engaged  first,  and  it  was  not  until  the  blade  load  reached  a 
value  of  26  pounds  that  the  upper  lobes  began  to  assume  load 
(Figure  54) . 
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Figure  52.  Photoelastic  Test  Model  of  Fir-Tree  Attachment 


Figure  53.  Photoelastic  Test  Setup  of 
Fir-Tree  Attachment  Model. 
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Figure  54.  Photoelastic  Fringe  Pattern  Illustrating 
Higher  Load  Carried  by  Lower  Lobes 
(Dark  Field)  . 
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When  both  sets  of  lobes  were  engaged,  the  fraction  of  the 
total  blade  load  carried  by  each  set  of  lobes  was  unknown. 
Hence,  an  arbitrary  average  stress  was  defined  across  the 
minimum  section  between  the  two  pairs  of  lobes  and  between 
the  upper  lobes  and  the  blade  proper  as  the  ratio  of  total 
blade  load  to  the  corresponding  minimum  section  area. 
Reference  is  made  to  the  stress  concentration  factor  (SCF) , 
which  is  defined  as  the  ratio  of  the  maximum  to  the  average 
stress,  as  an  "apparent"  or  "effective"  SCF.  Maximum 
stresses  in  all  cases  were  obtained  from  the  photoelastic 
patterns  at  each  load. 

Results  obtained  for  the  first  set  of  tests  for  the  blade 
section  between  the  upper  and  the  lower  lobes  are  given  in 
Table  XIII. 


TABLE  XIII.  SCF 

BETWEEN  UPPER  AND 

LOWER  LOBES 

Fringe 

Order 

Blade  Load 
(lb) 

Maximum 

Stress 

(psi) 

Average 

Stress 

(psi) 

Apparent 

SCF 

1/2 

16.3 

170 

75.8 

2.24 

1 

34.1 

340 

158.8 

2.14 

3/2 

58.4 

510 

272.0 

1.87 

2 

82.0 

68C 

382.0 

1.80 

Maximum  deviations  of  individual  load  measurements  from  the 
average  for  each  set  of  readings  at  a  particular  maximum 
stress  level  were  less  than  12  percent,  and  average  devia¬ 
tions  were  of  the  order  of  5  percent.  Results  given  in 
Table  XIII  apply  only  to  the  case  where  the  upper  lobes  do 
not  make  contact  with  their  mating  surfaces  until  the  blade 

load  exceeds  F  =  26.1  pounds  in  the  model, 
o  c 

Analogous  results  of  the  second  tests  for  the  section  between 
the  upper  lobes  and  the  blade  are  given  in  Table  XIV.  Load 
values  are  accurate  to  within  -rQ.5  pound,  and  fringe  values 
are  accurate  to  within  0.04  of  a  fringe. 


115 


TABLE  XIV.  SCF 

BETWEEN  UPPER 

LOBES  AND 

BLADES 

Fringe 

Order 

Blade  Load 
(lb) 

Maximum 

Stress 

(psi) 

Average 

Stress 

(psi) 

Apparent 

SCF 

1.00 

38.6 

340 

114 

3.  CO 

1.37 

56.4 

‘  466 

166 

2.81 

1.88 

♦ 

o 

CO 

640 

236 

2.71 

l . — : - - — - - -  - -  - - j 

Typical  fringe  patterns  from  which  the  data  were  obtained 
are  shown  in  Figures  55  and  56-  Again,  it  should  be  noted 
that  the  results  given  above  are  applicable  only  to  tolerances 
existing  in  the  dimensions  of  the  fabricated  model  used. 


Figure  55.  Typical  Dark-Field  Photoelastic  Pattern 
at  High  Load. 
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% 

Figure  56.  Light-Field  Photoelastic  Pattern 
Corresponding  to  Internal  Load 
Range. 


When  the  full  blade  load  was  carried  by  the  lower  lobes  for 
loads  less  than  or  equal  to  F  ,  and  further  loads  were  shared 
by  the  upper  and  lower  lobes, °then 


SCF  -  K  = 

Cl 


a  A 
in 


+'C(P  “ 


V 


for  F  >  F. 


(32) 
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and 


K  =  -g-  for  F  <  F  (33) 

S  F  O 


where  F  =  total  blade  load 

A  =  cross-sectional  area  (0.215  sq  in.) 

C  =  fractional  load  carried  by  lower  lobes  when 
both  lobes  were  active  at  all  times 

a  =  maximum  stress  at  section  for  blade  load,  F 
m 

For  fringe  order  1/2,  F  <  Fq  and  Equation  (33)  were  appli¬ 
cable,  and  for  this  case  the  apparent  and  true  stress- 
concentration  factors  were  identical.  For  the  remaining 
cases.  Equation  (32)  may  be  written  in  the  forms 


K  F  +  C(F_  -  F  )  =  (a  A\  (3! 

a  l  o  2  o  j  (  *2  J 

where  the  subscripts  1,  3/2,  and  2  refer  to  the  data  for  the 
corresponding  fringe  order.  By  eliminating  C  from  the  above 
equations,  three  equations  for  K  were  obtained: 

a 


f~  <F2  -  Fx)  (37) 
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r~ 


m. 


/  T~» 

U'  /  *\ 


li.- 


1/2  iFl  -  r°’] 


A_ 

r. 


u3/2  “  Fl) 


(38) 


A 

F 


(F2  "  F3/25 

(39) 


which  yielded,  respectively,  K  •-  2.34,  2.34,  and  2.28.  The 

a 

average  value  obtained  from  the, four  sets  of  data  gave 
K  =  2.30.  By  using  this  value  of  K  (2.30)  in  Equations  (34), 

S  cl 

(35),  and  (36),  C  -  0.714,  0.670,  and  0.674,  respectively, 
for  an  average  value  of  C  --  0.686. 


A  fairly  comprehensive  picture  of  the  stress  levels  expected 
to  develop  at  the  blade  section  between  the  upper  and  lower 
lobes  for  a  variety  of  fits  at  the  fir-tree  joint  has  evolved. 
Thus,  Equation  (32)  may  be  written 


C  A 

y  _ Q!_ 

eff  '  F 


(40) 


where  R  =  Fq/F  =  fraction  of  total  load  at  which  the  upper 
lobe  first  became  engaged. 

The  effective  stress  concentration  factor  as  defined  by 
Equation  (40)  varied  linearly  with  R  from  a  value  of  1.59 
at  R  =  0  to  2.30  at  R  =  1.0.  The  maximum  stress  developed 
in  the  prototype  at  operating  speed  (i.e.,  blade  load  of 
2500  pounds)  is  given  as  a  function  of  R  in  Figure  57. 

The  fractional  load  between  appearances  of  fringe  orders  1, 

3/2,  and  2  was  divided  between  the  upper  and  lower  lobes 

in  the  ratio  C:(l-C).  Thus,  there  was  an  independent  measure 

of  K  ,,  at  R  =  0.  The  averaae  value  so  determined  was  K  = 
etr  err 

1.52  as  compared  with  the  1.59  value  calculated  previously. 

A  similar  calculation  for  K  the  effective  stress  concen- 

ef  f 

tration  factor  between  the  upper  lobes  and  blade,  gave  an 
average  value  for  R  =  0: 


2.46 


With  simultaneous  engagement  of  both  sets  of  lobes,  a  maximum 
stress  of  22,600  psi  was  expected  at  this  section  in  the 
prototype  (beryllide)  blade. 

The  results  of  this  photoelastic  stress  analysis  of  the 
fir-tree  connection  between  the  blade  and  hub  of  the 
2400 °F  turbine  wheel  indicated  the  following: 

1.  With  simultaneous  engagement  of  all  lobes  of  the 
fir  tree,  the  lower  lobes  would  carry  68  percent 

of  the  total  blade  load.  This  fraction  of  the  total 
load  may  be  reduced  slightly  in  the  prototype  because 
of  the  camber  on  the  outer  boundary  of  the  outer  hub 
lobes. 

2.  The  maximum  stress  between  the  upper  and  the  lower 
blade  lobes  (23,100  psi)  was  about  2  percent  higher 
than  that  between  the  upper  lobe  and  blade  proper 
(22,600  psi)  at  R  =  0.  This  difference  was  less 
than  the  possible  experimental  error.  By  using 
the  experimental  value  Kef^  =  1.52  rather  than 

the  experimental  value  1.59,  a  value  2  percent 
lower  (rather  than  2  percent  higher)  was  obtained. 
With  the  change  in  stiffness  due  to  the  camber 
mentioned  above,  it  was  probable  that  for  R  =  0 
the  maximum  stress  would  be  slightly  higher  between 
the  upper  lobe  and  the  blade  than  elsewhere  in  the 
fir-tree  joint. 

3.  The  results  formed  a  basis  for  a  comprehensive 
experimental  stress-analysis  program  leading  to 
improved  design  and  establishment  of  allowable 
tolerance  limits  for  inspection. 

4.  It  further  appeared  to  be  desirable  to  investigate 
stress  conditions  resulting  from  nonsymnetric  en¬ 
gagement  of  upper  and  lower  lobes.  Ibis  could  also 
be  accomplished  through  the  proper  use  of  shims. 

In  addition,  analysis  of  the  data  indicated  that  the  first 
blade  model  either  behaved  nonlineariy  at  low  loads  or  had 
initial  edge  stresses  introduced  during  machining.  Further, 
the  most  recent  m.easurgjsQnt?  of  the  modulus  of  elasticity 
of  beryllide  had  lowered  the  model  load  for  similitude  with 
the  prototype  to  this  range,  which  made  behavior  at  low  loads 
of  even  more  significance.  Therefore,  a  new  blade  model  was 
made,  and  the  behavior  under  low  loads  was  investigated. 
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The  stress  distribution  in  the  elements  of  the  fir-tree  con¬ 
figuration  under  load  depended  markedly  on  the  deformation 
accompanying  the  load.  The  maximum  stresses  in  the  blade 
varied  with  the  load  in  a  nonlinear  and  discontinuous  manner 
in  the  model  tests. 

Steel  shim  stock  was  placed  between  the  disk  and  the  upper 
blade  lqbes  to  ensure  loading  of  only  the  upper  lobes.  The 
load  was  slowly  applied,  and  the  fringe  order  at  the  fillet 
above  the  upper  lobes  was  determined  at  approximately  every 
5-pound  increment  of  load,  to  a  maximum  total  of  35  pounds. 

During  loading,  it  was  noticed  that  the  increment  of  fringe 
order  accompanying  the  load  increase  from  22  to  27  pounds 
was  of  a  larger  magnitude  than  the  fringe  order  over  equal 
load  increments.  Therefore,  the  load  was  reduced  from  35 
to  approximately  22  pounds  to  give  an  additional  check 
measurement  of  the  fringe  order  at  that  load.  All  of  the 
load  was  then  removed,  and  the  specimen  fringe  order  was 
again  zero,  as  at  the  start  of  the  load  cycle. 

The  test  data  were  then  reduced  to  give  absolute  fringe  order 
as  a  function  of  load,  and  the  results  were  plotted  as  shown 
in  Figure  58.  The  data  defined  two  apparently  separate 
curves.  Equally  surprising  was  the  difference  in  fringe 
order  at  a  22-  to  23 -pound  load  obtained  on  the  ioading- 
versus-unloading  cycle. 

The  only  explanation  for  the  phenomena  described  above  is  that 
the  load  distribution  at  the  loaded  blade-disk  interface 
changes  suddenly  (i.e.,  discontinuously}  at  a  load  somewhere 
between  22  and  27  pounds  on  loading  and  at  some  load  below 
23  pounds  on  unloading.  The  discontinuity  is  caused  by 
sticking  and  slipping  of  one  contacting  surface  against  the 
other. 

There  was  again  some  evidence  of  the  stick-slip  process  when 
the  shims  were  moved  to  give  contact  at  the  lower  blade  lobes 
only.  However,  this  occurred  at  smaller  load  increments,  and 
it  was  possible  to  plot  a  continuous  curve  for  maximum  stress 
versus  load  as  shown  in  Figures  59  and  60,  These  latter 
curves  also  indicated  bending  of  the  blade  in  addition  to  the 
tension.  This  is  shown  by  the  difference  in  stress  at  the 
two  fillets  between  the  lower  and  upper  blade  lobes. 
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FRINGE  ORDER,  1(1.0  I  -  340  PSI) 
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Figure  58.  Model  Stress  Above  Upper  Blade  Lobes 
(Load  on  Upper) . 
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Figure  59. 
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Figure  60.  Model  Stress  Between  Upper  and 

Lower  Blade  Lobes  {Load  on  Lower) . 
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The  nonlinearity  of  maximum  fillet  stress  as  a  function  of 
load  is  due  to  the  influence  of  elastic  deformations  on  stress 
distribution  in  the  fir-tree  configuration.  Thus,  in  addition 
to  geometrical  similitude  of  the  photoelastic  model  and  its 
prototype,  there  must  be  similitude  in  load  distribution  if 
photoe.iastic  data  are  to  be  applicable  to  the  turbine  itself. 


When  the  relationship  between  the  applied  model  load  and 

the  prototype  load  was  derived,  it  was  recognized  that, 

should  deflections  influence  load  distribution,  a  further 
problem  could  arise  with  regard  to  its  applicability  for 
the  following  reason:  in  the  model,  both  the  disk  and  the 
blade  were  of  the  same  material  (CR-39) ,  having  a  value 
E  •-  3  x  106  psi;  while  in  the  prototype,  the  disk  was  of 

IN- 100  (E  =  24  x  10s  psi)  and  the  blade  was  of  beryllide 
P 

(Ea  =  39  x  106  psi).  This  became  a  real  problem  when  inter- 
P 

preting  the  results  obtained  in  terms  of  prototype  stress 
levels  to  be  expected. 


In  contrast,  the  result  for  the  disk  (Figure  61)  is  singularly 
linear.  The  probable  explanation  for  this  is  thSlT  during  the 
slipping  process,  all  of  the  relative  motion  occurred  between 
the  blade  and  the  shim,  so  that  the  load  distribution  on  the 
disk  model  remained  ostensibly  unchanged. 


The  redistribution  of  load  in  the  prototype  may  be  expected 
to  differ  from  that  in  the  model  for  two  reasons:  (1)  the 
difference  in  elastic  modulus '-of  the  blade  and  disk  materials 
and  (2)  the  vibration  of  the  engine  during  loading,  which  may 
eliminate  the  stick-slip  process  (the  change  in  load  distri¬ 
bution  would  be  a  continuous  process) . 


Despite  the  limitations  of  the  model  studies,  useful  design 
information  may  nevertheless  be  obtained  from  the  photoelastic 
results. 


The  prototype  stresses  were  calculated  at  the  three  critical 
fillet  sections  in  the  blade  and  disk  on  the  assumption  that 
both  the  blade  and  disk  material  was  IN-100  and  recalculated 
on  the  assumption  that  the  material  was  beryllide.  These 
stresses  were  plotted  as  functions  of  the  percent  of  full 
load  (2500  pounds)  carried  by  the  upper  and  lower  lobe  pairs 
in  Figures  62  and  63.  The  discontinuity  shown  in  Figure  58-. 
was  "faired"  out  in  developing  the  upper  curves  in  Figure  631 
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Figure  63.  Maximum  Blade  Stress  Betv.een 
Upper  and  Lower  Lobes  and 
Above  Upper  Lobes. 
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The  curves  of  Figure  63  are  labeled  with  regard  to  similitude 
with  a  material  having  E  =  39  x  10e  psi  (beryllide)  or  one 
having  E  -  24  x  106  psi  (IN- 100).  For  the  case  of  a  beryl¬ 
lide  blade  and  an  IN-100  disk,  the  blade  stresses  would  be 
higher  than  those  shown  by  the  lower  curves,  and  use  of  the 
upper  curves  for  design  purposes  would  not  be  overly  conserva¬ 
tive. 

PHOTOELASTIC  TEST  CONCLUSIONS 

Review  of  the  photoelastic  test  results  revealed  that  the 
indicated  stresses  were  somewhat  lower  than  those  calculated 
during  the  fir-tree  attachment  design  studies.  However,  as 
reported  in  the  design  studies  (Paragraph  3. 3. 1.2. 2),  there 
were  a  number  of  minor  differences  between  the  fir-tree  at¬ 
tachment  modeled  in  the  photoelastic  tests  and  the  fir-tree 
attachment  analyzed  in  the  design  studies.  Thus,  care  was 
exercised  when  comparing  the  stresses  from  the  two  studies. 

The  major  contribution  of  the  photoelastic  test  was  the  data 
on  the  relative  effects  of  tolerances  and  load  sharing  between 
the  lobes.  This  information  made  the  photoelastic  testing  a 
valuable  tool  in  establishing  the  final  fir-tree  configu¬ 
ration. 


Figure  64.  Loading  Considered  in  Fir-Tree  Analysis. 
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APPENDIX  IV 

THERMAL-CONDUCTION  TEST  PROGRAM 


GENERAL 


The  tests  described  in  this  appendix  were  conducted  to 
determine  the  thermal  conductance  across  the  interface  of 
the  beryllide  and  IN-100  contacting  surfaces  in  the  two- lobe 
fir-tree  disk/blade  attachment. 

CONDUCTION  RIG  DESIGN 


The  instrumented  test  specimens  consisted  of  tv/o  IN-100 
cylinders  0.75  inch  in  diameter  over  a  length  of  1.5  to  2.0 
inches  (Figure  65 ) .  The  upper  specimen  had  an  enlarged 
section  1.25  inches  in  diameter  at  its  upper  end,  which  was 
heated  by  an  induction  coil.  The  lower  specimen  terminated 
in  a  block  approximately  4  x  4  x  2.5  inches,  in  which  pro¬ 
vision  was  made  for  circulating  cooling  water. 

The  ends  of  the  specimens  away  from  the  enlarged  sections 
were  ground  to  a  finish  of  rms  20-24.  Thermocouples  were 
attached  to  the  specimens  along  the  heat  flew  path.  The  four 
thermocouples  on  the  circumference  and  the  thermocouple  on  the 
specimen  axis  at  axial  stations  0.125  and  0.875  inch  from  the 
interfaces  (i.e.,  the  specimen  ends)  served  to  give  a  measure 
of  the  uniformity  of  temperature  and  load,  over  the  specimen 
cross  sections.  Additional  thermocouples  0.5  inch  from  the 
interfaces  and  at  the  specimen  axis  were  used  in  conjunction 
with  the  others  to  measure  the  axial  temperature  distribution. 
The  thermocouple  installation  is  shown  in  Figure  66. 

With  these  two  specimens,  it  was  possible  to  directly  measure 
the  interface  resistance  to  heat  flow  or  its  reciprocal,  the 
thermal  conductance  of  the  interface,  as  a  function  of  the 
pressure.  However,  specimens  of  the  beryllides  could  not  be 
satisfactorily  instrumented  and  were,  therefore,  fabricated 
as  thin  disks  0.75  inch  in  diameter  and  0.0625  and  0.125  inch 
thick,  ttiese  disks  ware  ground-finished  on  both  sides  to  16 
rms  and  presented  two  interfaces  with  the  IN-100  specimens. 

For  comparison,  an  IN-100  disk  0.125  inch  thick  was  also 
fabricated  and  finished  to  the  same  rms  value. 

During  the  tests,  thermal  insulation  was  packed  around  the 
specimens  over  the  full  length  of  the  0.75-inch-diaraeter  sec¬ 
tions  while  the  specimens  were  in  alignment  and  under  a  com¬ 
pressive  load  of  about  200  pounds. 
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Figure  G5.  Conduction  Rig  Test  Setup. 
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CONDUCTION  RIG  TESTING 


Loads  were  applied  to  the  specimens  through  a  spherical  seat 
loading  head  by  a  hydraulic  test  machine,  with  power  to  the 
induction  heating  coil  supplied  by  an  oscillator.  Tests  ware 
conducted  at  interface  pressures  of  100,  1000,  and  2000  psi. 
Thermocouple  readings  were  taken  every  15  minutes  until 
steady- state  conditions  were  reached  for  each  load. 

The  thermal  conductivity  variation  with  temperature  for  IN-100 
is  given  in  Figure  6?  and  was  used  to  determine  the  heat  flow 
across  the  interface.  Similar  data  for  beryllide  are  also 
given  in  Figure  67  and  were  used  to  calculate  the  temperature 
drop  across  the  beryllide  disk. 

Temperatures  measured  on  the  circumference  of  the  specimen  at 
Station  1  were  40°  to»50°F  higher  than  those  on  the  axis.  A 
similar  difference  was  observed  at  Station  6,  but  with  tem¬ 
peratures  on  the  circumference ' lower  than  those  on  the  axis. 
Tills  was  due  to  the  shortness  of  the  transition  distances 
between  the  heat  source  and  heat  sink  and  the  instrumented 
test  length  of  the  specimen.  However,  the  variation  in  tem¬ 
peratures  as  measured  at  Stations  3  and  4  was  of  the  order  of 
10 °F  in  each  case,  and  the  average  temperatures  were  within 
a  degree  or  two  of  the  axis  temperatures. 

The  average  interface  temperature  in  the  tests  with  the  beryl¬ 
lide  disks  was  1109°F;  with  an  IN- 100  disk,  it  was  1064°F; 
and  with  IN-100  versus  IN-100  (a  single  interface) ,  it  was 
873 °F. 

Calculated  contact  resistances  based  on  the  test  data  are  in 
reasonable  agreement  with  the  sparse  published  values  avail¬ 
able  in  the  literature  for  interface  temperatures  and  pres¬ 
sures  considered  in  Table  XV. 

Differences  in  contact  resistance  among  the  various  test  con¬ 
figurations  were  not  .considered  to  be  large,  in  view  of . the 
individual  differences  between  the  characteristics  of  the  con¬ 
tacting  surfaces  in  each  test,  the  differences  in  relative 
orientation  of  tool  marks  on  contacting  surfaces,  etc. 

The  thermal  contact  resistance  between  IN-100  and  beryllide 
as  measured  in  these  tests  was  roughly  one-half  as  large  as 
the  value  assumed  in  the  design  phases  of  the  uncooled 
turbine. 
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TABLE  XV.  COMPARISON  OF  THERMAL  CONTACT  RESISTANCE  AS 
MEASURED  IN  REPORTED  TESTS  WITH  PUBLISHED 
DATA  FOR  SURFACES  OF  APPROXIMATELY  SIMILAR 
RMS  VALUES 


Contact  Resistance 


Interface  <Hr-FtV,F/Btu) 


Material 

RMS 

xn<-t=x  xa tre 

Temperature 

•F 

1000  psi 

2000  psi 

Steel/ steel* 

3-32 

180 

0.00045 

0.00040 

Steel/steel* 

3-16 

600 

0.00010 

0.00005 

Beryllide/j.N-100  favg) 

16-24 

1100 

0.00026 

0.00021 

IN-i00/lN-100 

16-24 

875 

0.00027 

0.00023 

*Kouwenhoven,  W.  R..  and  J.  H.  Potter,  THERMAL  RESISTANCE  OF 
METAL  CONTACTS,  The  Welding  Journal,  Vol.  27,  No.  10, 
October  1948,  ppV  515s-520s. 


CONDUCTION  TEST  CONCLUSIONS 


The  test  results  indicated  that  the  contact  resistance  was 
approximately  50  percent  smaller  than  that  estimated  and  used 
in  the  rotor  thermal  analysis  of  the  final  design  mechanical 
studies.  However,  the  estimated  contact  resistance  used  for 
the  thermal  analysis  was  so  small  that  further  reduction 
would  have  had  little  or  no  effect  on  the  calculated  disk 
temperatures . 


136 


Sggtrfiy  Ogwiftfattcw 


DOCUMENT  CONTROL  OaTA  -RAD 

'&*<  urltf  oi  tit f*.  kcfy  ct  +*4  mntroflic*  ctutt  fur  w-H+n  th*  *m*i'*/J 


ORIGINATING  ACtIVITV  t+  «vdt 9ij 

AiResearch  Manufacturing  Company  of  Arizona 
402  S.  36th  Street 

Phoenix.  Arizona  _ _ .  . _ _ _ 


J.  MCPCMT  TlTUC 

2400*F  UNC00LED  TURBINE/MATERIAL  PROGRAM 
VOLUME  II  -  AERODYNAMIC  AND  MECHANICAL  DESIGN 


4.  oeiCROTivI  Notts  (Tjrf  •ittpmt  m.4  Afei; 

Final  Renort 


Scnteaber  1970 


M  COMTHACT  OH  «MHT  HO. 

DA  44-177-AMC-183(T) 

A.  riojCC'  NO. 

Task  1GI62204A01409 


to  cttfmturio*>  «T«tCH<MT 


[M.  ORitlMATOirt  RKAOftT  KUtltXRtl) 


USAAVLABS  Technical  Report  70-32B 


This  document  is  subject  to  special  export  controls,  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  U.  S.  Army 
Aviation  Materiel  Laboratories,  Fort  Eustis.  Virginia  23604. 


•  t.  IRONI9NIN4  MLttASV  ACTIVITY 

U.S.  Army  Aviation  Materiel  Laboratories 
Fort  fcustis,  Virginia 
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Volume  II  of  a  3-voliunc  report 


If.  ABSTRACT 

This" report  describes  serodynamic/and  mechanical  design  studies  conducted  to  establish 
the  preliminary  design  of  a  2400* F  uncoolcd  turbine  component.  The  design  activity  was 
conducted  concurrently  with  materials  research  to  define  the  properties  of  IN-100  and 
AiRcsist  13  and  to  develop  a  high-temperature  turbine  material  from  intermetal lie 
compositions  of  beryllide  (Voli»t~l).  The  turbine  mechanical  design  studies  were  based 
on  estimates  of  IN-100  and  beryllide  material  properties. 

Preliminary  design  studies  were  conducted  to  define  a  number  of  aerodynamic  designs  for 
efficiency  comparisons  and  blade  stress  analyses,  leading  to  selection  of  a  candidate 
turbine  design.  Preliminary  thermal  and  centrifugal  stress  analyses  were  conducted 
for  each  disk  of  the  candidate  turbine. 

Tne  design  established  for  the  uncooled  turbine  was  a  three-stage  axial  flow  turbine 
with  equal  work  split  among  stages.  The  first-  and  second-stage  blades  and  vanes  were 
of  beryllide  material,  and  the  third-stage  rotor  was  to  be  integrally  cast  IN-100. 
Extensive  thermal  and  stress  analyses  were  conducted  on  first-stage  turbine  components 
to  evaluate  the  feasibility  of  the  uncooled  turbine  design. 

The  design  studies  were  terminated  when  it  became  apparent  that  the  proposed  beryllide 
turbine  material  lacked  sufficient  stress-rupture  strength  for  rotor  application  and 
did  not  possess  adequate  ductility  for  either  stator  or  rotor  application. 
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